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ABSTRACT 


Observations  of  a  highly  structured  discharge 
and  current  instability  in  unseeded,  shock  heated 
argon  flowing  through  a  region  of  large  transverse 
magnetic  field  are  presented  and  analyzed.  The  exper¬ 
imental  results  are  related  to  an  extension  of  the 
theory  of  ionization  (electrothermal)  instability  in 
partially  ionized  gases  in  which  the  presence  of 
multiply  charged  ions  is  taken  into  account.  A  critical 
magnetic  field  for  instability  onset  is  predicted. 

The  growth  rate  and  the  angular  dependence  of  the 
growth  rate  are  calculated  and  compared  with  experi¬ 
mental  observations.  An  enhanced  interaction  between 
the  current  discharge  and  the  plasma  flow  is  observed. 
The  magnitude  of  the  interaction  is  compared  with  the 
magnitude  of  the  Hall  parameter  in  the  plasma. 

A  rotation  of  the  current  filaments  to  form  a  helical 
structure  is  observed.  The  formation  of  this  structure 
is  predicted  by  the  use  of  a  simple  minimum  energy 


consideration . 
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CHAPTER  1  INTRODUCTION 


The  subject  of  this  work  is  the  investigation  of 
a  magnetically  induced  instability  in  a  partially 
ionized  plasma.  The  instability  manifests  itself  as 
a  highly  non-uniform  current  distribution  composed  of 
a  number  of  current  pinches  or  streamers  when  the  plasma 
is  subjected  to  crossed  electric  and  magnetic  fields 
(see  Figure  1.1). 

The  behavior  of  the  instability  is  explained  by  an 
extension  of  the  theory  of  the  ionization  or  electro¬ 
thermal  instability  in  which  the  presence  of  multiply 
charged  ions  is  accounted  for. 


1.1  A  Brief  Review  of  the  Literature 

The  discovery  of  the  ionization  instability 
mechanism  grew  out  of  the  work4 6 ' 47 ' 51 ' 9 1  directed 
toward  utilizing  extrathermal  or  nonequilibrium 
ionization,  produced  by  electric  fields  or  other  means. 


' 


-2- 


to  enhance  the  conductivity  of  the  working  fluid  of  an 
MHD  generator.  The  instability  mechanism  was  first 
proposed  by  Velikhov  and  Dvkhne^^  ana  Kerrebrock^ . 

Both  authors  showed  that  the  system  of  equations  consist¬ 
ing  of  a  generalized  Ohm's  Law,  an  electron  energy 
balance  equation  for  Te>Tp,  and  an  ionization  equil¬ 
ibrium  relation  (Saha's  equation),  is  unstable  to 
small  perturbations  in  electron  density  and  temper¬ 
ature.  Velikhov  and  Dykhne  derived  a  relation  for  a 
critical  Hall  parameter  above  which  instability  occurs 
and  showed  that  the  plasma  effective  conductivity  will 
be  appreciably  lowered  when  the  instability  is  present. 
Kerrebrock,  restricting  his  analysis  to  the  case  of 
slight  ionization  and  cons  cant  electron  collision 
frequency,  derived  the  growth  rate,  frequency,  and 
wavelength  of  the  oscillations  and  showed  that  the 
growth  rate  exhibits  a  directional  dependence. 

The  first  experimental  investigations  into  the 
nature  of  the  instability  were  performed  by  Nedospasov 
and  Shipuk^  and  Klepeis  and  Rosa^ .  Nedospasov  and 
Shipuk  used  a  pulsed  discharge  in  a  mixture  of  argon  and 
mercury  vapor  placed  between  the  poles  of  an  electro¬ 
magnet.  They  found  that  the  Hall  field  exhibited  a 
saturation  and  ceased  to  depend  on  magnetic  field. 
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and  that  the  time  averaged  conductivity  decreased 
inversely  with  Hall  parameter.  Klepeis  and  Rosa, 
using  a  disc  geometry  MHD  generator  duct  with  argon 
seeded  with  cesium,  also  observed  a  decreased  Hall 
voltage  as  well  as  the  formation  of  highly  luminous 
radial  spokes  at  large  values  of  magnetic  field. 

This  investigation  was  continued  and  extended 
by  a  number  of  workers7 ' 45 ' 101 ' 114 ' for  example. 

It  was  shown  that  striations  or  streamers  form  in  the 
plasma  when  the  critical  Hall  parameter  is  exceeded. 

The  streamers  are  inclined  at  a  characteristic  angle 
to  the  main  current. 

The  theory  of  the  instability  was  considered 
further  by  Solbes^-88,  Zampaglione^-^ ,  Dykhne^,  Velikhov 
et  al-^5,  Nelson7*^  '  7  ^  '  7  8  ,  Hiramoto8^,  and  Gilpin28. 

All  the  authors  considered  perturbation  solutions 

“!  I  C 

to  the  above  mentioned  equations.  Velikhov  et  alXL-> 
and  Nelson7 ^  used  numerical  methods  to  investigate 
the  effect  of  the  boundary  conditions  due  to  the  gener¬ 
ator  duct  walls  and  electrodes  on  the  behavior  of 
the  instability.  Gilpin28  related  the  instability 
onset  conditions  to  the  breakdown  condition  for  a 
gas  discharge. 
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FIGURE  1.1(a) 


Flow 

Direction 


Photograph  of  discharge 
without  transverse  magnetic  field 


(  B=0,  Ey=40  v./cm. ;  produced  by  applied  potential; 
1-6000  amps,  (frame  5  ) ;  10  ysec.  interframe  delay; 
exposure  times  per  frame  are  0.5,  0.5,  0.5,  0.2, 

0.1  ysec.  respectively.) 


Increasing 
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Flow 

Direction 


Figure  1.1(b) 

Photograph  of  Discharge 

with  Transverse  Magnetic  Field 
2 

(B  --  1.63  w/m  ;  produced  by  motion  through 

magnetic  field;  -  50  v/cra;  I  -  1000  amps; 

2  ysec  inter frame  delay;  0.2  ysec  exposure  time 
per  frame) 


' 


-6- 


The  instability  behavior  in  a  flow  of  gas  was 
investigated  by  Gilpin  and  Zukoski27 ' 126 ,  Kerrebrock  and 
Dethlef sen52 ,  Riedmiiller87  ,  Vitshas,  Golubev  and  Malikov117, 
Shipuk  and  Pashkin102,  Golubev  et  al  31 ,  and  Malikov66. 

Gilpin  and  Zukoski,  Riedmiiller,  Vitshas  et  al, 
and  Shipuk  and  Pashkin  observed  regular  streamers 
oriented  at  an  angle  to  the  main  current.  The  effects 
of  Hall  voltage  saturation  and  increase  in  turbulence 
and  amplitude  of  fluctuations  with  increase  in  magnetic 
field  were  also  observed. 

O  1 

Golubev  et  al  and  Malikov00,  using  a  sodium  seeded 
argon  flow  observed  a  number  of  plasma  streamers  forming 
in  the  presence  of  a  magnetic  field.  Because  the  streamers 
were  oriented  primarily  along  the  Faraday  or  v*S  direction 
rather  than  at  an  angle  to  the  Faraday  direction,  the 
authors  were  reluctant  to  attribute  the  filamentary 
structure  to  the  presence  of  the  ionization  instability, 
despite  commonly  observed  behavior  such  as  the  increase 
in  fluctuation  amplitude  with  Hall  parameter  and  satur¬ 
ation  of  the  Hall  field  with  increasing  Hall  parameter. 
However,  it  is  shown  in  this  work  (Chapter  3)  that, 
during  the  formation  period  of  the  current  streamers, 
they  are  inclined  at  the  characteristic  angle  to  the 
^*6  direction  but  that,  due  to  constraints  on  the  current 
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f low  (constrained  to  flow  primarily  in  the  vxB  direction) , 
the  fully  developed  filaments  must  be  aligned  primarily 

-y 

in  the  v*B  direction. 

Zauderer121 ' 124  has  observed  current  pinches  in  a 
linear  MHD  channel  which  persisted  in  the  absence  of 
a  magnetic  field.  No  increase  in  turbulence  and  no 
saturation  of  Hall  voltage  were  observed.  It  thus  appears 
that  the  current  pinching  was  due  to  the  thermal  pinch 
mechanism  that  occurs  in  a  classical  high  pressure  arc, 
and  not  to  the  magnetically  induced  instability  invest¬ 
igated  here. 

Klingenberg"  observed  current  pinches  in  shock 
heated  gas  in  the  presence  of  a  magnetic  field.  He 
attributed  their  formation  to  electrode  effects  and  did 
not  consider  the  role  of  the  magnetic  field  in  producing 
a  current  instability.  The  value  of  Hall  parameter  in 
his  experiment  was  less  than  or  just  equal  to  the  minimum 
value  for  appearance  of  the  instability.  Under  these 
conditions,  small  variations  in  Joule  heating  rate  intro¬ 
duced  by  varying  the  electrode  arrangement  may  serve  to 
suppress  the  instability.  Another  factor  is  the  role 
played  by  impurities  in  his  experiment.  His  test  section 
was  constructed  of  plexiglass  and  the  interaction  of 
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the  hot  gas  with  this  plastic  surface  almost  certainly 
contaminated  the  shock  heated  gas.  It  has  been  recog¬ 
nized^  that  a  few  tenths  of  a  per  cent  of  molecular 
impurity  can  absorb  a  large  proportion  of  the  electron 
energy  and  thus  lower  the  electron  temperature.  This 
lowering  of  the  electron  temperature  will  damp  out  the 
instability.  In  the  present  work,  it  has  been  observed 
that  a  high  impurity  level  caused  by  the  failure  to  flush 
the  shock  tube  with  clean  gas  before  each  test  is  suffi¬ 
cient  to  produce  an  overall  brightening  of  the  gas, 
masking  the  instability. 


1.2  The  Nature  of  the  Present  Work 

All  of  the  previous  work,  in  which  bhe  authors 
recognized  the  existence  of  the  magnetically  induced 
instability,  has  been  conducted  in  low  temperature, 
seeded,  atomic  gases,  in  which  the  dominant  particle 
interaction  mechanism  was  electron  atom  collisions 
and  the  electron  density  was  almost  completely  due  to 
the  ionization  of  the  seed  atoms.  The  purpose  of  the 
present  work  is  to  consider  the  behavior  of  the  instability 
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in  the  much  more  general  and  widely  applicable  case  of 
an  unseeded,  partially  ionized.  Coulomb  collision  domin¬ 
ated  plasma  in  which  the  electron  density  is  due  directly 
to  the  ionization  of  the  atomic  gas. 

The  partially  ionized  gas  is  produced  by  shock 
heating  argon  as  described  in  Chapter  2 .  The  magnetic 
field  is  perpendicular  to  the  flow  direction  of  the  gas 
and  the  electric  field  is  produced  by  the  motional  inter¬ 
action  between  the  gas  flow  and  the  magnetic  field. 

In  Chapter  3,  the  theory  of  the  ionization  instab¬ 
ility  is  extended  to  the  Coulomb  collision  dominated  case. 
The  presence  of  multiply  charged  ions  is  accounted  for 
by  the  use  of  an  approximate  method  of  calculation  due 
to  Raizer^’ .  Expressions  for  the  dependence  of  the  onset 
of  the  instability  on  Hall  parameter  and  directional 
dependence  of  the  growth  rate  on  Hall  parameter  are 
derived.  -Numerical  calculations  are  made,  for  the  case 
of  shock  heated  argon,  of  the  temperature  nonequipartition , 
the  instability  onset  and  growth  rate,  and  the  angular 
dependence  of  the  growth  rate.  The  observed  experimental 
behavior  of  the  streamers  shows  good  agreement  with  the 
calculations.  The  breakup  of  the  streamers  into  smaller, 
secondary  filaments  is  related  to  the  angular  dependence 


of  the  growth  rate. 
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In  Chapter  4,  the  behavior  of  the  fully  developed 
current  streamers  is  examined.  An  enhanced  interaction 
between  the  plasma  streamer  and  the  flow  of  shocked 
argon  is  observed,  based  on  a  consideration  of  the  force 
balance  on  the  streamer.  It  is  suggested  that  this  large 
interaction  is  due  to  turbulence  produced  in  both  the 
streamer  and  argon  flow  by  the  instability  mechanism  of 
Chapter  3.  This  theory  is  supported  by  an  observed  corre¬ 
lation  between  the  degree  of  interaction  and  the  Hall 
parameter . 

A  rotation  of  the  plasma  streamers  to  form  a  heli¬ 
cal  structure  is  observed.  The  formation  of  this  structure 
is  explained  by  the  existence  of  a  component  of  magnetic 
field  parallel  to  the  axis  of  the  helix.  The  angle  of 
pitch  and  the  diameter  of  the  helix  are  predicted  on  the 
basis  of  a  simple  minimum  energy  consideration.  The 
predicted,  values  show  good  agreement  with  experiment. 

In  Chapter  5,  suggestions  are  given  for  possible 
applications  of  the  present  work.  Possible  future 
areas  of  research  are  indicated. 
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CHAPTER  2  DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS 


AND  DIAGNOSTIC  TECHNIQUES 


2.1  Method  of  Production  and  Properties  of  the 

Shock  Heated  Gas 

The  shock  wave  was  produced  by  a  pressure  driven, 
diaphragm  shock  tube  with  helium  heated  by  an  electric 
discharge  as  the  driver  gas.  The  electric  discharge  was 
produced  by  a  20,000  volt,  8700  Joule  capacitor  bank. 

The  driver  section  consisted  of  a  5  cm.  length  of 
5  cm.  diameter  by  1  cm.  thick  stainless  steel  tube. 

A  stainless  steel  flange,  "o"-rings,  and  mechanical 
clamps  were  attached  to  one  end  of  the  tube  to  accomodate 

_  o 

a  mylar  polyester  film  (1.5  *  10  inch  thickness) 
diaphragm.  A  flange,  "o"-rings,  and  a  1  inch  thick 
plexiglass  insulator  were  fastened  to  the  other  end  of 
the  tube.  The  discharge  occurred  between  a  brass 
electrode  passing  through  the  center  of  the  insulator 
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and  the  inner  wall  of  the  stainless  steel  tube.  The 
driver  was  flushed  with  helium  and  filled  to  5  p.s.i. 
above  atmospheric  pressure  before  each  test. 

The  shock  tube  consisted  of  a  5  cm.  diameter  stain¬ 
less  steel  tube  in  which  the  shock  wave  was  allowed  to 
flow  for  a  distance  of  about  3  meters  before  reaching 
the  magnetic  interaction  region.  The  unshocked  gas  was 
at  a  pressure  of  0.5  torr  and  room  temperature.  Shocks 
of  up  to  Mach  15  with  a  sample  of  shock  heated  gas  5  to 
6  cm.  long  were  produced  3  meters  from  the  driver. 

The  equilibrium  properties  of  the  shock  heated  gas 
were  calculated  by  de  Leeuw^.  Appendix  1  gives  proper¬ 
ties  of  shock  heated  argon,  derived  from  de  Leeuw, 
for  the  pressure  and  Mach  number  range  used  in  this 
experiment.  A  shock  tube  of  similar  design  operated 
under  similar  experimental  conditions  was  studied  by 
Phillips8'1'.  He  experimentally  verified  that  the  proper¬ 
ties  of  the  shock  heated  gas  are  well  described  by  the 
usual  Rankine-Hugoniot  shock  jump  equations.  He  also 
showed  that,  at  the  initial  pressures  and  Mach  numbers 
used  here,  the  shock  heated  sample  had  uniform  properties 
over  the  5  to  6  cm.  sample  length,  within  the  limits  of 
his  diagnostic  methods.  Phillips  also  showed  that  the 
thickness  of  the  ionization  relaxation  region  is  much  less 
than  the  total  thickness  of  the  shock  heated  gas. 
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2.2  Description  of  the  Magnetic  Interaction  Region 

The  shock  tube  in  the  interaction  region  consisted 
of  a  60  cm.  length  of  5  cm.  diameter  pyrex  tubing,  beside 
which  two  large  magnetic  field  coils  were  placed.  The 
transverse  magnetic  field  was  produced  by  discharging  a 
5000  volt,  21,000  Joule  capacitor  bank  through  the  two 
coils,  series  connected  and  arranged  in  the  Helmholtz  pair 
configuration.  The  period  of  the  discharge  was  about  800 
ysec.  ,  much  longer  than  the  duration  of  passage  of  the 
shock  heated  gas.  The  coils  consisted  of  0.5  mm.  by 

1.5  cm.  brass  strip  spirally  wound  and  embedded  in  epoxy 
resin.  A  5  cm.  diameter  hole  was  provided  in  the  center 
of  the  coils  to  permit  viewing  of  the  passage  of  the 
shock  heated  gas.  Magnetic  fields  of  up  to  2  w./m.2 
were  produced  over  the  5  cm.  electrode  region.  The  field 
was  spatially  uniform  to  within  15  %  over  the  area  of 
the  viewing  port  and  was  constant  in  time  to  within  10% 
for  the  duration  of  passage  of  the  shock  heated  gas. 

Figure  2.1  shows  details  of  the  experimental  arrange¬ 
ment  in  the  magnetic  interaction  region.  Brass  electrodes, 

2.5  cm.  in  diameter,  were  machined  and  mounted  flush 
with  the  interior  wall  of  the  pyrex  section  of  the  shock 
tube.  The  electrodes  were  connected  together  by  a  20  cm. 
length  of  0.15  mm.  by  1  cm.  brass  strip. 
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2.3  Diagnostic  Methods 

The  main  diagnostic  tool  used  was  an  electronically 
shuttered,  image  converter  camera  (TRW  Model  ID),  per¬ 
mitting  observations  of  the  temporal  and  spatial  devel¬ 
opment  of  the  discharge  structure.  Using  the  camera, 
measurements  were  made  of  the  growth  rate  of  the  streamers, 
the  direction  of  maximum  growth  rate,  and  the  effect 
of  increasing  magnetic  field  strength  on  the  streamer 
structure.  In  addition,  the  velocities  of  the  streamers 
and  the  luminosity  front  and  the  dimensions  of  the 
streamers  were  measured.  Piezoelectric  pressure  trans¬ 
ducers  were  also  used  to  measure  shock  speed. 

The  total  electrode  current  was  determined  by  meas¬ 
uring  the  voltage  drop  across  a  section  of  the  brass 
strip  connecting  the  two  electrodes,  the  resistance  of 
which  was  accurately  known. 

Much  electronic  equipment  for  the  amplification  and 
delay  of  high  voltage  pulses  was  developed  in  order  to 
trigger  the  capacitor  banks,  oscilloscopes,  and  electronic 
camera  at  the  appropriate  times.  This  equipment  w ill  not 
be  described  here  since  it  is  of  reasonably  standard 
design.  However,  one  technique  developed  here  is  of  broad 
interest  to  many  workers  in  experimental  investigations 
of  plasmas  and  will  be  described  below. 
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One  problem  frequently  encountered  in  plasma 
research  is  due  to  the  large  amount  of  electromagnetic 
radiation  produced  whenever  a  high  energy  discharge  occurs, 
for  instance,  whenever  a  high  voltage  arc  is  struck  or 
whenever  an  energy  storage  capacitor  bank  is  discharged. 
This  radiation  invariably  activates  the  low  level  trig¬ 
ger  circuitry  on  any  nearby  oscilloscopes.  Frequently, 
such  as  in  the  operation  of  an  electrically  driven 
shock  tube  or  an  electrically  excited  laser,  the  operation 
of  the  oscilloscope  must  be  delayed  for  a  short  period 
of  time  in  order  to  observe  the  phenomena  of  interest. 
Although  one  delayed  time  base  can  be  obtained  with  a 
dual-time  base  oscilloscope,  in  many  cases  it  is  necess¬ 
ary  to  have  two  or  more  delayed  time  bases  to  observe 
phenomena  of  short  duration.  This  .is  impossible  to 
accomplish  unless  some  method  is  found  to  prevent  the 
unwanted  triggering  of  the  primary  time  base. 

One  modification  to  accomplish  this  purpose  will  now 
be  described.  It  is  directly  applicable  to  a  Tektronix 
555  oscilloscope  but  the  same  principle  may  be  used  on 
any  oscilloscope  of  similar  design. 

The  time  base  generators  incorporate  a  "lockout 
multivibrator"  (Figure  2.2)  which  prevents  operation 
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of  the  time  base  under  certain  conditions;  it  is  the 
circuit  which  permits  single  sweep  operation  and  the 
delayed  sweep  of  the  secondary  time  base.  Figure  2.2 
shows  a  circuit  modification  designed  to  prevent  the 
operation  of  the  primary  time  base  (time  base  "A") 
until  an  external  high  level  pulse  (-30  volts)  is 
applied.  This  external  pulse  either  causes  the  time 
base  to  sweep  immediately  or  removes  the  lock-out  and 
permits  the  next  low  level  trigger  signal  to  trigger 
the  time  base. 

When  switch  Sg  is  closed,  the  time  base  triggers 
immediately  upon  application  of  the  -30  volt  pulse 
(note  that  the  M sweep  function"  switch  must  be  set  at 
"single  sweep";  if  it  is  in  any  other  position,  the  sweep 
free-runs  and  cannot  be  synchronized  or  triggered) . 

When  Sg  is  open,  the  time  base  is  enabled  (the 
ready  light  comes  on)  by  the  -30  volt  pulse.  The  time 
base  may  then  be  triggered  by  the  usual  means.  When 
Sg  is  open,  the  time  base  also  operates  normally  in  every 
position  of  the  "sweep  function"  switch,  provided  no 
-30  volt  pulses  are  applied. 

The  delayed  -30  volt  pulse  can  be  obtained,  for 
instance,  from  a  thyratron  delay  unit,  which  can  be 
made  relatively  insensitive  to  electromagnetic  inter¬ 
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CHAPTER  3  ONSET  OF  THE  IONIZATION  INSTABILITY 


3.1  Outline  of  Theory 

In  a  partially  ionized  gas  subjected  to  an  electric 
field,  the  possibility  exists  of  a  non-equipartition  of 
energy  between  electrons  and  heavy  particles.  Specifi¬ 
cally,  we'  would  expect  an  excess  of  the  electronic  ther¬ 
mal  energy  over  the  heavy  particle  thermal  energy, 
where  the  electronic  thermal  energy  is  essentially  deter¬ 
mined  by  gain  in  the  electric  field  and  loss  due  to 
collisional  processes,  and  the  electron  concentration  is 
determined  by  approximate  ionization  equilibrium  at  the 
electron  temperature.  As  in  the  work  of  Zauderer-*-^ , 
Kerrebrock^  '  ^  ,  and  Zukoski  and  Gilpin^  r 1^6  ,  we  repy 
on  convection  to  maintain  the  gas  temperature  nearly  cons¬ 
tant  while  the  gas  flows  through  a  test  region. 


.  . 
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Thus,  we  consider  a  situation  where  the  electrons 
comprise  a  gas  that  is  weakly  coupled  thermally  to  the 
heavy  particle  gas,  but  very  strongly  coupled  to  the 
valence  electrons  of  the  gas,  and  exhibiting  a  Maxwellian 
distribution  of  energies.  One  more  parameter,  then,  is 
required,  in  addition  to  the  thermodynamic  variables  of 
the  gas,  in  order  to  determine  the  state  of  the  plasma. 

The  electron  temperature  is,  of  course,  the  fundamental 
variable  that  must  be  added.  In  the  following  sections, 
a  method  of  calculation  of  the  electron  temperature  is 
given . 

The  question  of  the  stability  of  this  non-equilibrium 
state  will  be  considered  in  detail.  In  particular,  the 
theories  of  Velikhov  and  Dykhne41^,  and  Kerrebrock49 
indicate  that  a  unique  mechanism  of  ionization  instability 
exists  in  a  non-equilibrium  plasma  placed  in  crossed 
electric  and  magnetic  fields.  This  instability  manifests 
itself  as  a  constriction  of  the  discharge  into  current 
streamers.  The  theory  is  extended  to  the  high  temperature, 
multiply  ionized.  Coulomb  collision  dominated  regime. 

The  predictions  of  the  extended  theory  will  be  compared 
with  the  experimental  results  obtained  in  partially 
ionized  argon  subject  to  crossed  electric  and  magnetic 


fields . 
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3.1.1  Calculation  of  the  Electron  Temperature 

The  rate  at  which  energy  is  gained  by  the  electrons 
is  the  scalar  product  of  the  electron  current  ^  and  the 


electric  field  in  the  frame  of  the  moving  gas,  E. 
Ohm's  law,  neglecting  diffusion  and  electron  inertia 


.  90 

effects,  is 


(3.1) 


ctE  -  0)Tj  x 


3 


B 


where  a  is  the  electron  scalar  conductivity, 
w  is  the  electron  cyclotron  frequency, 
and  i  is  the  electron  mean  free  time  or  the 
reciprocal  of  the  electron  collision  frequency. 

In  a  coordinate  system  as  shown  in  Figure  3.1, 
Ohm's  Law  is: 
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or,  for  j  =  0,  (no  Hall  current  flows). 


°0Ex  =  "Uy 


■’y  = 


°oEy 


Thus , 


Figure  3.1 
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For  a  Maxwellian  electron  distribution  function,  the 

energy  loss  of  the  electrons  per  elastic  collision  with 

heavy  particles  is  2me/mK  times  the  difference  in  energy 

of  electrons  and  heavy  particles^ ' ^ ^  ,  where 

is  the  electron  mass  and  mK  is  the  mass  of  the  K'th 

species  of  heavy  particles.  To  account  for  possible 

inelastic  collisions,  a  correction  factor  <$  is  intro- 

K 

duced  which  is  close  to  unity  for  monatomic  gases  (elastic 

collisions)  but  which  is  10  to  10  for  diatomic  and 

polyatomic  molecules  in  which  vibrational  and  rotational 

3  9  9  2 

modes  may  be  excited  '  '.  The  energy  balance  of  the 

electrons  is  then3^ ' ' ^ ^ ^  ^ ^ 

°0Ey  =  E (neveK>  <2me/mK)6K(3k/2) (Te-TK)  ...(3.3) 
K 

where  n  is  the  number  density  of  electrons,  veK  is  the 
electron  collision  frequency  with  the  K'th  species  of 
heavy  particles,  k  is  Boltzmann's  constant,  Te  is  the 
electron  temperature,  TR  is  the  temperature  of  the  K'th 
species.  The  TR's  are  generally  identical  and  are  equal 
to  the  gas  temperature  Tg  ^ . 
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In  the  plasma  of  interest  here,  (partially  ionized 
argon),  there  are  two  species  of  heavy  particles,  argon 
atoms  and  argon  ions.  Thus,  the  energy  balance  equation 
becomes : 


a 


ne (2me/ma) (3k/2) (Te-Tg) (vei  +  vea)  ...(3.4) 


since  6  is  equal  to  unity  for  both  species  of  heavy 
K 

.  .  49  92 

particles4 ' ' 


TT106 
Nov;  , 


ne  e‘ 


o 


. . .  (3.5) 


me<vei  +  vea> 


So,  the  energy  balance  equation  is: 


Ej  =  (3km|/mae2) (Te-Tg) (vei  +  Vea)  ...(3.6) 

Thus,  given  an  (electric  field  in  the  frame  of  the  gas) 
and  given  vei  and  vea,  we  can  calculate  the  electron 
temperature . 

In  the  temperature  regime  of  interest  here,  (Mach  8 

6  0 

to  Mach  16  argon  shocks) ,  vei  is  much  greater  than  vea  ; 
Coulomb  collisions  are  the  dominant  collision  process. 
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vei  is  given  by  Rose  and  Clarke90  as: 

Z  ne  InA 

vei  —  “  ”  "  P~  TTp  ...(3.7) 

1.51  x  106  T3/Z 

e 

Now,  we  also  need  a  relation  between  ne ,  the  electron 
density,  and  Te,  the  electron  temperature  or  mean  thermal 
energy.  This  relation  is  given  by  the  ionization  equil¬ 
ibrium  relation  or  Saha's  equation. 

3.1.2  Equilibrium  of  Electrons 

As  noted  in  the  introduction  to  this  chapter,  we 

may  think  of  the  free  and  valence  electrons  as  comprising 

a  single  gas.  Electron-electron  collisions  tend  to  force 

this  gas  toward  a  Boltzmann  distribution,  whereas  collisions 

of  electrons  with  atoms  or  ions  and  radiation,  in  general, 

4  9  51 

tend  to  perturb  it.  As  noted  by  Kerrebrock  '  and  Ben 
Daniel  and  Tamor9,  there  is  a  value  of  electron  density 
or  ionization  fraction  below  which  the  electrons  are 
influenced  more  by  energy  exchange  with  heavy  particles 
than  by  mutual  energy  exchange.  At  electron  densities 
sufficiently  far  above  this  critical  value,  Coulomb 
collisions  between  electrons  override  all  other  processes 
and  lead  to  a  Maxwellian  distribution. 
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Following  Kerrebrock49 ,  the  condition  for  a  Maxwellian 
distribution  of  electron  energies  is: 


>  > 


(2  me/ma) (nennoen  +  neni0ei) 


(3.8) 


where  aee,  aen'  aei  are  cross  sections  for 

electron-electron,  electron-neutral,  and  electron-ion 
collisions  respectively;  ne,  nn,  and  nj_  are  the  electron, 
neutral,  and  ion  number  densities  respectively;  me  and 
ma  are  the  electron  and  atom  (ion)  masses. 

Using  ne  =  Zna  and  dividing  by  nenncee,  we  have: 


Now ,  since  me  <<  ma  and  (ae^/aee)  -  1  (ref.  90), 
we  have : 


(ne/nn)  >>  2  (me/ma) (aen/oee) 
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Now,  aen  =  10  20  m2  (ref.  20),  a0e  =  5  x  10  17  m2  (ref.  90), 
and  ine/ma  -  1.36  x  10  .  Using  these  numerical  values, 

we  expect  the  electrons  to  have  a  Maxwellian  distribution 
at  the  electron  temperature  whenever 

(ne/nn)  >>  5  x  10“9  ...(3.9) 

Thus,  for  the  range  of  temperatures  and  densities 
encountered  in  this  experiment,  we  would  expect  the 
electrons  to  have  a  Maxwellian  distribution  of  energies 
about  the  electron  temperature. 

It  remains  to  consider  the  effects  of  radiation  from 
the  electron  gas  on  the  electron  energy  balance. 

3.1.3.  Effect  of  Radiation  Losses  on  the  Electron 
Energy  Balance  and  the  Electron  Temperature 

The  purpose  of  this  section  is  to  compare  the  rate 
of  electron  energy  loss  due  to  radiation  with  the  rate  of 
energy  transfer  to  the  heavy  particles  due  to  electron- 
heavy  particle  collisions. 
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Petschek  et  al3®  show  that  the  main  energy  loss 


mechanism  for  shock  heated  argon  in  the  temperature  regime 
of  interest  here  is  continuum  radiation  due  to  free- 
bound  electron  transitions.  By  measuring  the  variation 
in  light  intensity  versus  distance  behind  the  shock,  the 
authors  were  able  to  show  that  this  cooling  of  the  gas 
was  closely  accounted  for  by  the  theory  of  continuum 

radiation  from  a  partially  ionized  gas  formulated  by 
Unsold3-^.  More  recently,  Horn,  Wong,  and  Ber shader 3^ 
considered  the  contribution  of  emission  from  spectral 
lines  to  the  total  radiative  emission  energy  loss.  They 
found  that  the  4p4s  and  3d4p  Ar  lines  contributed  rather 
large  losses,  with  the  result  that  the  total  loss  is 
approximately  equal  for  both  continuum  emission  and  line 
radiation . 

In  this  section,  we  will  consider  only  the  continuum 
radiation.  From  ref.  38,  the  total  radiative  energy  loss 
per  unit  volume  per  unit  time  for  optically  thin  continuum 
radiation  is: 


ergs/cm  3/sec 


. . . (3 . 10a) 


for  ne  in  cm 


3  and  T 


e 


0T.r 

m  K . 
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. . .  (3.10b) 

for  n  in  m“3,  T  in  °R. 
e  e 

The  collisional  loss  rate  is  given  by: 


pc  =  2(me/ma)  ne  <3k/2) (Te  -  Tg)  vei 


Using  eg. (3.7)  for 


=  3.73x10  3^  Z  In  A  (l - £-  watts/m3 

c  \  T  /  T1/2 

e  e 

...(3.11) 

Now,  the  ratio  of  the  collisional  loss  rate  to  the 
radiative  loss  rate  is: 


P 


c 


Q 


5.44x10 


16 


vg  +  (kTe/h) 


Z 

,T 

'eff 


InA 


. . .  (3.12) 


o  o  y 

Following  Horn  et  al  ,  we  take  Z^ff=  1.5,  and 

hv  =  2.85  eV.  In  addition,  we  take  Z  =  1 ,  InA  =  6, 

y 

kT  -  1  eV,  and  (1  -  (Tn/Tp) )  “  o.5,  as  typical  values 
e  y  c 

for  the  experimental  conditions. 

We  then  have : 


P  /Q  -  1 . 2 x 102 


. . .  (3.13) 
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We  see  that  the  collisional  loss  rate  is  about  two 
orders  of  magnitude  larger  than  the  radiation  energy  loss 
rate,  and  we  are  justified  in  neglecting  the  effect  of 
radiative  energy  loss  on  the  electron  temperature. 
Including  the  contribution  of  line  emission  to  the  total 
radiative  loss  as  described  by  Horn  et  al  does  not 
invalidate  this  conclusion. 


3.1.4.  Calculation  of  the  Electron  Density 


The  calculation  of  the  electron  density,  assuming 
ionization  equilibrium  at  the  electron  temperature, 
proceeds  from  Saha's  equation: 


. . .  (3.14) 


I  is  the  ionization  potential,  g+/ga  is  the  ratio  of 

statistical  weights,  TQis  the  electron  temperature, 

n  ,  n- ,  and  n  are  the  electron,  ion,  and  neutral  number 
e  i  a 

densities,  respectively. 
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As  we  shall  see,  under  the  conditions  of  this  exper¬ 
iment,  electron  temperatures  will  be  high  enough  to  ionize 
the  second  and  perhaps  third  levels  of  the  argon  ions. 

Under  these  conditions,  we  must  use  the  system  of  Saha 
12  5 

equations  : 


n  .  .  n^ 
m+1  e 


u. 


=  A 


E+i.^3/2 


n 


m 


u 


exp (  -Im+i/kT  )  ...(3.15) 


m 


where  nm+]_/  ne/  and  nm  are  the  m+1  multiply  ionized  ion, 

electron,  and  m  multiply  ionized  ion  number  densities 

respectively.  um+i/um  statistical  weight  ratio 

for  the  electronic  states  of  the  ions.  I  , .  is  the 

m+1 

ionization  potential  for  the  m'th  ion. 

One  can  readily  appreciate  that  the  calculations  for 
electron  density  based  on  this  non-linear  system  of 
equations  are  extremely  involved  and  time  consuming. 

For  each  pair  of  temperature  and  density  values,  a  non¬ 
linear  system  of  algebraic  equations  must  be  solved  to 
determine  the  concentrations  of  ions  of  different  charge 
and  hence  the  electron  concentration.  In  addition,  the 


incorporation  of  the  full  system  of  Saha  equations  in  the 
calculation  of  the  dispersion  relation  for  ionization 


. 
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instability  (see  following  sections)  would  cause  the 
resulting  expression  to  become  very  involved  and  greatly 
complicate  its  derivation. 

For  these  reasons,  we  will  use  an  approximate  method 

•  8  o 

of  calculation  of  the  electron  density  due  to  Raizer  and 

described  in  detail  in  Appendix  2.  This  approximation 

rests  on  two  assumptions.  The  first  is  the  assumption 

that  the  ion  number  density  nm  and  the  ionization  potentials 

I  +1  are  considered  to  be  continuous  functions  of  the  ionic 

charge  multiplicity  m,  obtained  by  connecting  the  discrete 

values  of  nm  and  I  by  continuous  curves.  We  denote  these 

two  continuous  functions  by  n (m)  and  I (m)  respectively.  The 

second  assumption  is  the  approximation  that  the  average 

value  of  the  ionic  charge  multiplicity,  also  the  average 

number  of  free  electrons  per  original  atom 


m  = 


j m  n (m)  dm  ne 

-  --  - -  .in  ■  ■  ■■  ■  1 1  i  i  .him  ~'ii  ~~~  «■ 

fn  (m)  dm  n 


. . . (3.16) 


is  exactly  equal  to  that  value  of  m  for  which  the  ion 
distribution  function  n (m)  has  a  maximum.  (In  the  above 
equation,  m  is  the  average  value  of  m,  n0  is  the  electron 
density,  ng  is  the  total  heavy  particle  density) . 


. 


-33- 


As  shown  in  Appendix  2 ,  the  above  assumptions  lead 
to  the  following  simple  transcendental  equation  in  ffl: 


or, 


where 


A  T3/2 
e 


exp 


I  (m) 


I (m+%) 


A  = 


3/2 


. . . (3.17) 


. . .  (3.18) 


Because  the  right  hand  side  is  a  logarithmic  function 
of  m,  two  or  three  succesive  approximations  are  sufficient 
to  obtain  a  fairly  accurate  value  of  m,  given  the  function 
I (m) .  As  shown  in  Appendix  2 ,  the  function  I (m)  may  be 
approximated  quite  accurately  for  the  first  four  or  five 
levels  of  argon  by  the  expression 


I (m)  =  15  m  electron-volts 


. . .  (3.19) 


3.1.5  System  of  Equations  for  Ionization  Instability 

In  this  section,  we  will  present  the  equations  nec¬ 
essary  to  derive  the  dispersion  relation  for  ionization 
instability.  Three  equations  are  necessary,  a  generalized 
Ohm's  Law  or  equation  of  motion,  an  energy  conservation 


■ 


: 
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or  energy  balance  law  for  the  electrons,  and  a  relation¬ 
ship  between  n  and  the  mean  electron  thermal  energy  (T  ) , 
i.e.,  Raizer's  approximation  to  the  set  of  Saha  equations. 

i)  Generalized  Ohm's  Law: 

Spitzer^*^  derives  the  following  generalized  Ohm's 
Law  from  the  equations  of  motion  for  ions  of  charge  +  Z  and 
electrons : 


mi  me 

Z  p  e2 

3 1 

1  / 

-I-  .  [ 

m  • 

eZp  ' 

i 

—  E  +  v  x 


-  Zme 


n: 


-  (m-  -  Zm„)  j  x  b 


*) 


.  .  .  (3.20) 


where  p  =  nimi  +  neme  and  n  is  the  electrical  resistivity. 
The  (3j/3t)  term  is  negligible  for  times  greater  than 

a  small  fraction  of  a  microsecond,  that  is,  inertial  effects 

,  .  •  i  -]  105,94 

are  negligible 

Mow,  for  n  =  (mev/nee^)  the  term 

(mj_  -  Zme)  (j  x  g)  /eZp  becomes  n(w/v)j  x  S/B,  where  w  is 
the  electron  cyclotron  frequency.  So,  Ohm's  Law  becomes: 


■ 


. 


' 


. 
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E'  +  vxB  +  —  (mi  VPe  -  Zme  VpJ 
eZp  \  / 

=  n  ^ j  +  jx (w/v) (B/B)^ 

Now,  P0  =  nekTe  and  P±  =  n^kT.;.  Thus,  VPp  =  k  (n^ VT0  +  T  Vn  ^ 
and  VPi  =  k  (neVTj  +  TiVne)  .  So,  the  term  (m.j_VPe  -  ZmeVPj,)/eZp 
becomes,  using  the  fact  that  me<<mi,  (k/e) (vTe  +  Te(vne/ne)). 
Thus,  Ohm's  Law  becomes: 


-x  wB 

=  j  +  j  x  _  ...  (3.21) 

vB 


where  a  =  1/n  =  n^e^/m^v,  and  6=2'  +  vxB  is  the 
0  c 

electric  field  seen  in  the  frame  of  the  moving  gas. 


ii)  Energy  balance  law  for  electrons: 


The  electronic  energy  balance  is  obtained  as  a 
moment  of  the  Boltzmann  equation^ ^ ^ .  The  result  is: 
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d 


dt 


e 

e 


+ 


£ 


j -E  -  neAee 


. . .  (3.22) 


where  ee  is  the  mean  electron  thermal  energy  (3kTQ/2) . 

£j_  is  the  ionization  energy  I  =  I  (m+^)  .  j’-S  is  the  energy 
input  rate  due  to  Joule  heating.  neAee  is  the  energy 
loss  rate  from  the  electrons  and  is  equal  to 
2  (me/m  )  nev  ( ee  -  £j_)  +  neQr  (see  Sect»  3.1.1).  Qr  is  the 
radiation  loss  which  we  have  shown,  in  Section  3.1.3,  to 
be  negligible.  v  is  the  collision  frequency  for  electron- 
ion  collisions.  So,  the  energy  balance  equation  becomes: 


(3k/2) (d (neTe) /dt)  +  (ne/ng) (di/dm) (dne/dt) 

+  I(dne/dt)  =  j-E  -  2 (me/ma)nev (3k/2) (Te  -  Tg) 

...  (3.23) 


iii)  Raizer's  approximation:  relationship  between  ne 
and  the  mean  thermal  energy. 

Raizer's  approximation  gives  (Section  3.1.4): 


3/2 

ne  =  mng  =  A  Te  exp 


I  (m+J$) 


kT, 


. . . (3.24) 


We  express  this  relation  in  differential  form  as 


follows : 


. 
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dm  =  (A/ng)  exp(-I/kTp) 


(3Te 


dTe/2  +  IdTe/kTg 


-  T^dl/kj 


Dividing  by  the  original  expression  and  rearranging: 


. . .  (3.25) 


3.1.6  Growth  Rate  of  the  Ionization  Instability 

In  this  section,  we  calculate  the  growth  rate  and 
the  angular  dependence  of  the  growth  rate  of  the  ioniz¬ 
ation  instability.  We  assume  that  small  initial  perturb¬ 
ations  in  the  quantities  nQ ,  2,  j,  and  T0  are  described 
as  follows: 

ne  =  neo  +  nel  exP  (  *K*r  -f  iffit  )  ...(3.26a) 

2  =  2q  +  exp (  i2*r  +  iffit  ) 


. . . (3.26b) 


■- 


■ 

i 
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. . . (3 . 26c) 


. . . (3.26d) 


We  assume  that  there  is  no  component  of  j  or  E  or 
variation  of  properties  in  the  direction  parallel  to  S; 


i.e.,  we  assume  that  E,  j,  and  all  variations  are  co-planar 


and  that  B  is  perpendicular  to  this  plane.  We  neglect 
terms  higher  than  first  order.  We  take  the  electron 
collision  frequency,  v,  as  given  by  equation  (3.7). 

We  take  Z  as  the  average  ionic  charge  multiplicity,  m, 
defined  in  Section  3.1.4. 

Inserting  the  assumed  form  of  the  perturbation,  we 
have,  to  first  order: 


1 


t0  +  (3/2)x0(Tel/Te0) 


T 


V 


. . .  (3.27) 


v 


T 


. . .  (3.28) 


Now,  since  a 


2  3/2 

(n0e  /mev) ,  a  «  ,  and,  to  first 


order , 


a 


a 
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The  generalized  Ohm's  Law  (equation  3.21),  neglecting 
the  diffusion  terms  for  the  moment,  becomes 


(• 


a  E  =  j  +[a)Tj  x  (B/B) 


) 


To  zero' th  order, 


a0^0  =  ^0  x  (B/B) 


(• 


) 


. . . (3.29) 


To  first  order,  ana  subtracting  the  zero'th  order 


equation , 


->  -> 


-> 


an  E,  +  (3T_.,/2t_J  =  3,  +  ( wx  Q  j  ±><  (B/B) 


0  1  oo  el  ^  eO 


B/B)j 


( 


+  “To  <3Tel/2TeO)  "  (nel/ne0 


>)  (v<S/B)) 


. . .  (3.30) 


The  energy  balance  equation,  eq.  (3.23),  is 

( 3kT  /2 )  (dn  /dt)  +  (3n0k/2)  (dTQ/dt) 
e  ^ 

+  (n  /n  )  (dl/drn)  (dne/dt)  +  I  (dne/dt) 
e  9 

=  j-E  -  3k(me/ma) (n0/x) (T0  -  Tg) 


* 
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To  zero' th  order. 


VE0  -  3k(meA»a)  (ne0/'o>  (Te0  ~  V  •••(3.31) 

To  first  order,  and  subtracting  the  zero'th  order 
equation , 


iffi  |  (3k/2)nelTe0  +  (3k/2)ne0Tel  +  in0l 


+  (di/dm) (neQnel/ng)  | 


n 


-  3  k (me/ma)  {  --- 


eO 


7ng)  j 

->  -> 

=  3l'E0 

+  VE1 

3T  n 
el 

1  1 

f-3 

CD 

O 

1 

H3 

1 _ 1 

^  ne0 

2T  n- 
eO 

f  ^neOTel^To^  +  ^nel/To^TeO  Tg^| 


. . .  (3.32) 


Mote  that  including  the  diffusion  terms  in  Ohm's  Law 
would  add  imaginary  terms  to  Ohm's  Law  and  thus  on  the 
right  hand  side  of  equation  (3.32),  through  the  Joule 

~y  ~y. 

heating  terms  j^*Eq  and  jg*E^.  fkese  terms  would  not 
affect  the  growth  rate  but  would  introduce  a  finite  phase 
velocity  (see  references  34  and  72) . 

The  ionization  equilibrium  expression,  equation  (3.25), 


becomes : 


. 


. 
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n 


+  — 


eO 


dl  \  n 


el 


kT  n  n  dm  /  n  . 
eO  g  /  eO 


We  note  also,  from  equation  (3.25),  that 


dlnT 


1  + 


ne0  dI 


kTeO  ngdm. 


(Tel/Te0> 


dlnn 


3  I 
2  kT 


eO 


(n  n/n  ) 
el  eO 


.  .  .  (3.34) 


We  need  two  more  conditions  in  order  to  eliminate 
the  first-order  quantities  from  the  dispersion  relation, 
equation  (3.32).  We  shall  exclude  electromagnetic  wave 
propagation  from  consideration  and  take  the  applied  field, 
B,  as  being  constant  in  time.  Thus, 


curl  E  -  0 


. . .  (3.35) 


We  also  make  the  plasma  approximation  that  there  is  no 
net  charge  density  and  thus, 


div  J=0 


. . .  (3.36) 
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We  now  use  equations  (3.29),  (3.30),  (3.31),  (3.33), 

(3.34),  (3.35),  and  (3.36)  to  eliminate  the  first  order 

quantities  from  the  dispersion  relation  (see  Appendix  3 
for  a  detailed  description) .  The  expression  for  the 
growth  rate  as  given  in  Appendix  3  is  then: 

(3  dln  Te  ,  , 

y  =  i(3  =  —  (s.inzc|>  -  cosz<j>) 

a nn  ,.1  (2  din  n^ 
u  eO  e 


+ 


_2wTns.in<}>  coscj) 


3  din  Te  Te0  din  To  ) 

2  din  n0  Te0”Tg  dln  ne  ^ 


(  3  kTeO 

•IT" 


1  + 


dln 

1 

<d| 

^  ! 

+  Hsl 

1 

dl 

dln 

ne  " 

ng 

I 

din 

+  1 


-1 


. .  .  (3.37) 


where  ( ir /2 )  ~<f>  is  the  angle  between  the  perturbed  current, 
j ^ ,  and  the  unperturbed  current,  jQ.  (dln  Te/dln  ne)  is 
given  by  equation  (3.34)  .  I  =  I(m+^)  ,  where  in  is  the 
average  ionic  charge  from  equation  (3.18) .  (di/dm)  is 
given  by  : 


dl  dl (m  +  h) 
dm  dm 


d 

dm 


15m  +7.5 


15  eV. 


for  argon  (see  Section  3.1.4). 
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•  2  ■ 

Jq/oq  is  the  Joule  dissipation  in  the  gas  (see 
equation  3.2) .  u)t0  is  the  Hall  parameter  and  aQ  is  the 
scalar  conductivity  calculated  using  the  electron  temp¬ 
erature  T  and  the  average  ionic  charge  multiplicity  in. 

As  shown  in  Appendix  3 ,  the  angle  <j>  for  which  the 
growth  rate  is  a  maximum  is  given  by: 


r. 

'2 


tan 


-1 


3  din  T 

e 

2  din  ne 


.  .  .  (3.38) 


Using  this  value  of  <{>  in  the  growth  rate  expression, 
equation  (3.37),  one  obtains: 


. .  .  (3.39) 
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Wo  see  that  the  growth  rate  increases  with  increasing 

oTq  and  also  depends  on  the  relation  between  equilibrium 

electron  density  and  electron  temperature  through  the 

term  (din  T  ) / (din  n  ). 

e  0 

We  also  note  the  dependence  on  the  temperature 
difference  between  the  electrons  and  heavy  particles, 
through  the  term  Te0/  (T0Q-Tg)  .  As  TeQ  approaches  T  , 
the  growth  rate  becomes  strongly  negative,  and  any 
perturbation  will  be  quickly  damped. 

3.1.7  Physical  Mechanism  of  the  Ionization  Instability 

The  following  physical  explanation  of  the  ionization 

...  73 

instability  is  due  to  Nedospasov  .  The  physical  mechan¬ 
ism  can  be  understood  by  considering  the  fluctuation  of 
the  Joule  heating  connected  with  the  fluctuation  of  the 
plasma  density. 

Let  us  assume  that  in  the  shaded  region  of  the 
figure  (Figure  3.2),  the  electron  density  is  larger  than 
in  the  surrounding  space  by  a  small  amount  6ne.  We 
assume  that  the  main  current  is  directed  along  the 
y-axis  and  the  magnetic  field  is  perpendicular  to  the 
plane  of  the  figure  and  directed  out  of  the  page. 


' 
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Figure  3 . 2 


From  the  generalized  Ohm's  Law,  we  have: 


E 


E 


Ox 

Oy 


CO  T 


For  simplicity,  we  neglect  the  fluctuations  of  t.  The 

component  of  current  normal  to  the  streamer  is  conserved 

since  we  take  div  J  =  0.  Therefore,  electric  charges 

appear  on  the  boundary  of  the  streamer  and  change  the 

electric  field  inside  the  streamer.  We  note  that  the 

86 


85 


tangential  field  cannot  change 


since  we  take  curl  £  =  0 


■ 


, 


The  additional  power  released  in  a  unit  volume  of 


the  streamer  due  to  Joule  heating  is 


6W . .  = 

6  j 

E  + 

i  6E  / 

jh 

Jn 

on 

Joj.  ±  ' 

since  6En 

-* 

=  0  and 

63j. 

=  0. 

Now , 

°3„  =  3 

e/ne^ 

since  j„  «  nQ 

<3„  =  enevii 

and,  from 

Figure 

3.2, 

-•oil 

jQsin  <J>  .  So, 

5 j.i  =  bsin  * 

Now,  from  Figure  3.2, 

E  =  En  sin  <j>  +  ooxE  cos  <j)  , 

on  °Y  oy 

and  E  =  i ~/o  from  the  generalized  Ohm’s  Law  in  the 
oy  ° 

absence  of  Hall  currents.  So, 


From  Figure  3.2,  jQJL  -  jQ  cos  $  •  Also,  since 
E,  «  1/a  for  j.  constant  and  o  n e,  then  a  1/n^. 

JL 


- 


. 
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Thus  ,  6E± 


-EOJL  (  6ne/ne)  .  From  Figure  3.2, 


Eox  =  Eoy  cos  * 


WTEoy  sin  ^ 


and  EQy  =  j Q/ a  as  above.  So, 


and 


Thus , 


E 


o 


01 
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cos  <J>  -  uxsin  <j> 
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SE. 
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'joxsin  <f>  -  cos  cj) 


\  6ne 
/  n^ 
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jh 


6n^  \  j 


jQ  sin  cj> 


le  \  ->o 


ne  /  a 


+  |  jQ  COS  (j) 
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( 


sin  <J>  +  w t cos  <j) 


wTsin  <|)  -  cos  (J) 
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e 
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Simplifying , 


5wjh  = 


j2  , 

J°  ,;„2 


in^ 4>.  -  cos^l  +  2wxsin  <{>  cos  <J> 
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Now,  the  increase  in  the  energy  lost  due  to  collisions 
with  the  gas  atoms  is 


since  the  loss  rate  is  proportional  to  both  nQ  and 
Te~Tg  or  simply  TQ  if  is  large,  (again,  we  neglect 
the  fluctuations  of  t  with  6Te  and  6ne)  . 

Thus,  we  see  that  in  a  streamer  with  Sne  >  0,  the 
increase  in  the  energy  released  due  to  Joule  heating  can 
exceed  the  increase  in  the  energy  lost  to  collisions  with 
the  heavy  particles.  This  will  occur  for  values  of  wx 
above  a  certain  value.  The  condition  6Wj^  >  5WC  also 
depends  on  the  relation  between  one/ne  and  <5Te/Te. 

The  greater  the  quantity 

*Sne/ne  6  In  ne 

——  or  - . . .  , 

STe/Te  6  In  Te 

the  more  6W-h  exceeds  6WC.  The  angular  dependence  of  the 
instability  can  also  be  seen.  It  is  easily  shown  that 
SWj^  is  a  maximum  for  an  angle  <j>  of  htan  ^wx  which  ap¬ 
proaches  xt/4  for  idt  large. 


SLy  :.i 
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In  summary,  we  see  that  the  instability  occurs  via 
the  following  mechanism:  in  a  streamer  with  6ne  >  0, 
the  increase  in  the  electron  energy  due  to  Joule  heating 
can  exceed  the  increase  in  the  energy  lost  due  to  collisions 
with  the  heavy  particles.  This  in  turn  increases  Te  and 
the  degree  of  ionization,  leading  to  a  further  deviation 
of  the  electron  density  from  the  equilibrium  value. 

The  neglect  of  tiie  fluctuations  in  collision  frequency 
is  not,  strictly  speaking,  justifiable  in  the  Coulomb 
collision  dominated  regime.  However,  the  inclusion  of 
the  collision  frequency  fluctuations  leads  to  the  same 
qualitative  result. 


- 


. 


50- 


3.2  Application  of  Theory  to  Experiment 


3.2.1  Calculation  of  the  Plasma  Parameters  for  Partially 
Ionized,  Shock  Heated  Argon 


To  obtain  the  steady  state  plasma  parameters  (n  ,  T 
oot,  a,  (din  T0)/(dln  ne)  )  needed  to  calculate  the  growth 
rate  and  critical  angle,  <j>  ,  the  following  equations 

were  solved  numerically: 


->■ 


30  +  “to3ox(b/b) 


. . . (3.29) 


^  ^  nu  n  3k 

j  •  E  =  2— (T  -T  ) 

Jo  o  0  eo  g 

raa  T  0  2 


. . . (3.31) 


n^  =  m  n 
eo  g 


3/2 

A  T  exp 
eo 


I  [m+h) 

kT 


eo 


. . . (3.24) 


->  -> 

h-Eo 


°0Ey 


3o/°0 


. . . (3.2) 


If  one  neglects  electrode  loss  and  the  potential 
drop  across  the  (very  small)  external  resistance,  E^, 
the  electric  field  in  the  gas  (see  Figure  3.1)  should 
be  equal  to  ^2®*  However,  electrode  losses  can  be  apprec 
iable .  To  eliminate  this  inaccuracy,  the  measured  values 


. 
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of  Ey  were  used  in  the  calculations.  Specifically,  it 
was  found  that  the  measured  values  of  were  closely 
approximated  by  the  function 


E 

y 


U2B  exp 

u 


. . . (3.40) 


where  B  =  0.312  w/m2(see  Figure  3.3). 

The  initial  values  of  plasma  density,  n  ,  and  temp¬ 
erature,  Tg,  of  the  shock  heated  gas  were  taken  from 
59 

de  Leeuw  and  are  given  in  Table  3.1. 

Figure  3.3(a)  shows  calculated  curves  of  electron 

temperature  versus  magnetic  field  strength.  The  maximum 

in  Tq  and  the  subsequent  decrease  with  increasing  B  is  due 

to  the  decrease  in  the  measured  values  of  E  with  increas- 

Y 

ing  magnetic  field  (Figure  3.3).  This  has  the  effect  of 
introducing  a  maximum  in  the  Joule  heating  rate  and 
consequently  a  maximum  in  the  electron  temperature 
through  equation  (3.31). 

The  steady  state  temperature  was  not  measured  directly 
in  this  work,  partly  because  the  onset  of  the  instab¬ 
ility  is  so  rapid  as  to  render  very  difficult  any  obser¬ 
vation  of  the  steady  state  conditions.  However,  Vasil' eva 
et  ai111'112  have  measured  the  electron  temperature  in 
shock  heated  argon  moving  into  a  magnetic  field.  They 
used  a  Mach  8.5  shock  propagating  into  argon  at  10  mm  Hg 


■ 


. 


-52- 


initial  pressure  and  in  the  presence  of  a  magnetic  field 
of  from  3000  to  6000  gauss.  Due  to  the  higher  gas  pres¬ 
sures  used,  their  Hall  parameter  was  below  the  critical 
value  for  instability.  They  were  thus  able  to  observe 
a  steady  state  electron  temperature  and  their  measure¬ 
ments  are  shown  in  Figure  3.3(a).  Their  measured  values 
lie  below  the  calculated  curves;  however,  they  observe 
substantial  nonequilibrium.  The  reason  for  the  lower 
measured  values  is  the  larger  collision  frequency  (due  to 
higher  pressure)  under  the  conditions  of  their  experi¬ 
ment.  Thus,  the  transfer  of  energy  from  electrons  to 
ions  is  more  efficient  and,  for  a  constant  Joule  heating 
rate,  the  gas  cannot  support  as  large  an  energy  non- 
equipartition . 

3.2.2  Variation  of  Hall  Parameter  with  Applied  Magnetic 
Field 

Calculations  based  on  equations  (3.29),  (3.31),  (3.24), 

(3.2),  and  (3.40)  give  the  change  in  the  plasma  parameters 
with  increasing  magnetic  field.  Figure  3.4  gives  the 
variation  of  the  Hall  parameter,  wt,  with  magnetic  field. 

We  see  that  tox  increases  with  increasing  B.  Tnus ,  from 
equation  (3.39)  ,  the  growth  rate  of  the  ionization  inotab 
ility  will  increase  with  increasing  B. 


We  also  note  that 


■ . 
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Variation  of  E  vzith  Applied  Magnetic  Field 
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Table  3.1 

Plasma  Parameters  for  Shock  Heated  Argon 
(thermal  equilibrium  assumed; 
initial  Ar  pressure  is  0.5  torr) 


-*■ 


PARTICLE  TEMPERATURE 


Figure  3.3  (a) 


Electron  Temperature  vs. 
Magnetic  Field  Strength  for 
Various  Shock  Mach  Numbers 
( A  Measurements  of  Vasil' eva  et  al) 
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as  the  Mach  number  is  reduced,  oot  increases.  This  refl¬ 
ects  the  fact  that  at  lower  Mach  numbers,  the  temperature 
is  lower,  electron  density  is  lower,  and  t,  the  mean 
free  time  for  electron  collisions,  is  then  much  greater. 
Thus,  we  see  that  the  growth  rate  of  the  ionization  instab¬ 
ility  increases  as  the  Mach  number  is  reduced,  at  least 
in  the  range  of  Mach  numbers  considered  here,  for  a 
constant  initial  temperature  and  pressure  ahead  of  the 
shock. 

3.2.3  Variation  of  Grow7th  Rate  with  Applied  Magnetic 
Field  and  Hall  Parameter 

Equations  (3.29),  (3.31),  (3.24)  ,  (3.2),  (3. 40), and 

(3.39)  allow  one  to  calculate  the  growth  rate  of  the  ioniz¬ 
ation  instability  versus  magnetic  field  and  versus  Hall 
parameter.  Results  of  these  calculations  are  presented 
in  Figures  3.5  and  3.6. 

We  see  from  Figures  3.5  and  3.6  that  there  is  a 
critical  magnetic  field  and  a  critical  Hall  parameter  above 
which  the  growth  rate  for  the  instability  is  positive. 

The  critical  magnetic  field  is  between  0.4  and  1.4  webers/ 
m^  and  increases  with  the  Mach  number  of  the  argon  shock 
wave.  The  reason  for  the  increase  of  the  critical  field 
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with  Mach  number  is  that  at  the  higher  Mach  numbers,  the 
equilibrium  temperature  and  density  are  higher,  thus  t, 
the  mean  free  time  for  electron  collisions,  is  lower 
and  a  greater  value  of  magnetic  field  is  needed  to  pro¬ 
duce  the  same  value  of  ui,  the  Hall  parameter. 

From  Figure  3.6,  we  see  that  the  critical  Hall 
parameter  has  a  value  in  the  range  of  two  to  three  and 
is  only  slightly  dependent  on  shock  velocity. 

In  Figures  3.7  and  3.8,  the  e-folding  time  (the 
reciprocal  of  the  growth  rate)  is  plotted  versus  magnetic 
field  strength  for  argon  heated  by  a  Mach  12  shock  wave. 
Measurements  were  made  of  the  e-folding  time  from  image 
converter  photographs  made  as  described  in  Chapter  2. 

The  length  of  an  individual  streamer  was  assumed  to  be 
proportional  to  exp  (t/x)  where  x  is  the  e- folding  time. 
The  length  of  the  streamer  was  measured  at  two  different 
times  (corresponding  to  two  adjacent  frames  of  the  five- 
frame  image  converter  photographs) .  Then  the  e-folding 
time  is  given  by  the  expression 

fc2  “  fcl 

X  =  . . — - 

In  (L2/L1) 

where  L2  and  are  the  streamer  lengths  measured  at  the 
times  t2  and  t±  respectively.  Measurements  made  in  this 
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Growth  Rate  vs.  Hall  Parameter 
for  Various  Values  of  Mach  Number 
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Figure  3.7 

e-Folding  Time  vs.  Hall  Parameter 
for  a  Mach  12  Argon  Shock 
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Figure  3.8 

e-Folding  Time  vs.  Magnetic  Field  Strength 
for  a  Mach  12  Argon  Shock 
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manner  are  shown  on  Figures  3.7  and  3.8.  The  measured 
e-folding  time  is  well  within  an  order  of  magnitude  of  the 
predicted  e-folding  time.  This  agreement  is  somewhat 
fortuitous  in  the  sense  that  the  theory  predicts  a  growth 
rate  for  small,  plane  wave  perturbations  whereas  the 
experimental  measurements  are  made  on  a  large  amplitude, 
highly  non-linear  disturbance  wherein  the  growth  rate  is 
certainly  influenced  to  a  large  extent  by  such  things  as 
convective  heat  transfer  and  radiative  energy  transport. 
However,  one  would  still  expect  the  growth  rate  to  be 
of  the  same  order  as  that  predicted  by  the  theory  and 
this  expectation  is  borne  out  by  the  measurements. 

3.2.4  Variation  of  the  Critical  Magnetic  Field  and 

Critical  Hall  Parameter  with  Shock  Mach  Number 

The  values  of  magnetic  field  and  Hall  parameter  for 
which  the  growth  rate  of  the  ionization  instability  is 
zero  are  plotted  in  Figures  3.9  and  3.10.  In  addition, 
the  values  of  magnetic  field  and  Kali  parameter  which 
produce  a  positive  e-folding  time  equal  to  the  flow 
duration  time  of  the  shock  heated  gas  are  plotted. 

Data  obtained  from  the  experiments  with  shock 
heated  argon  as  described  in  Chapter  2  are  also  plotted. 
Figures  3.9  and  3.10  show  clearly  that  the  theory  coriect 
ly  predicts  the  onset  of  the  instability. 
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3.2.5  Angular  Dependence  of  the  Streamer  Formation 


As  shown  in  Appendix  3,  the  growth  rate  of  the  ioniz¬ 
ation  instability  has  an  angular  dependence  relative  to 
the  unperturbed  current.  The  direction  of  maximum  growth 
rate  is  at  an  angle 


h  tan 


-1 


03  T 


3  d  In  T 


2  d  In  n 


e  J 


. . . (A. 3 .36) 


where  the  value  of  the  tan  ^  function  is  taken  between  0 
and  tt/2.  a  is  the  angle  between  the  average  or  unperturbed 
current  and  the  instability  or  streamer  current.  In 
Appendix  3,  it  is  also  shown  that  the  streamer  may  be 
produced  on  either  side  of  the  main  current  j  at  the 
characteristic  angle  a  (see  Figures  A. 3.1  and  A. 3. 2). 

Thus,  the  theory  indicates  that: 

1)  the  streamers  should  be  oriented  at  an  angle 

a  to  the  main  current  (the  main  current  is  in 

~y  #  , 

the  Faraday  or  vxB  direction) . 

2  8 

2)  as  pointed  out  by  Gilpin  ,  the  streamers 

will  tend  to  break  up  at  higher  Kail  parameters, 

forming  smaller  streamers  inclined  at  an  angle 

a  to  the  original  streamer. 
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Using  the  experimental  arrangement  described  in 
Chapter  2,  photographs  were  made  of  the  discharge  structure 
in  the  center  of  the  magnetic  field  region  as  the  shock 
heated  gas  passed  through  the  magnetic  field.  Figure 
3.11(a),  a  sequence  of  photos  separated  in  time  by  2  ysec. , 
shows  the  structure  formation.  In  particular,  the  third 
frame  shows  a  number  of  parallel  streamers  oriented  at  an 
angle  of  about  40  degrees  to  the  vx§  direction.  This 
agrees  with  the  theoretically  calculated  behavior  (see 
Figures  3.13  and  3.14).  As  the  structure  develops,  one 
streamer  appears  to  predominate,  conducting  most  of  the 
current  in  one  channel  or  arc  (see  frames  4  and  5  of 
Figure  3.11(a)  ).  The  structure  of  the  streamers  then 

becomes  oriented  primarily  along  the  Faraday  or  v*i§  direc¬ 
tion,  (see  Figure  3.11(b)  ).  This  orientation  has  been 

observed  previously  '  by  workers  at  the  Kurchatov 
Institute  in  the  USSR  using  a  low  temperature,  seeded 
plasma.  This  change  in  orientation  as  the  structure 
develops  is  very  likely  due  to  the  imposition  of  bound¬ 
ary  conditions  on  the  plasma  in  the  experimental  configur¬ 
ation,  i.e.,  to  the  fact  that  the  streamer  current  is 
constrained  to  flow  primarily  in  the  v*)^  direction  because 
of  the  electrode  configuration.  This  fact  is  not  taken 
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into  account  in  the  theory  presented  in  Section  3.1,  since 
the  streamers  are  assumed  to  be  small  perturbations  to  the 
main,  uniform  vxB  current.  The  theory  will  cease  to  des¬ 
cribe  the  exact  physical  situation  when  the  perturbation 
or  streamer  current  becomes  large  compared  to  the  main 
current . 

Figure  3.11(b)  clearly  shows  the  breakup  of  the  current 
filaments  or  streamers,  forming  smaller  streamers  at  an 

p  o 

angle  to  the  original,  as  predicted  by  Gilpin  .  As  shown 
in  Figure  3.12,  the  smaller  streamers  form  at  an  angle  to 
the  main  streamer  of  about  40  degrees. 

The  angle  a,  between  the  perturbed  current  and  the 
original  current,  was  calculated  as  described  in  Section 
3.2.1.  Figure  3.13  gives  curves  of  the  angle  versus  mag¬ 
netic  field  for  various  shock  Mach  numbers.  Figure  3.14 
shows  data  taken  from  a  large  number  of  image  converter 
photographs  of  the  streamer  structure.  The  solid  curve  is 
the  theoretical  curve  for  initial  gas  conditions  produced 
by  a  Mach  12  shock  wave  propagating  into  argon  at  0.5 
torr  initial  pressure. 
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Figure  3.11(a) 

Image  Converter  Photographs 
of  the  Discharge  Structure 
in  the  Center  of  the  Magnetic  Field  Region 

Exposure  duration  =  0.2  ysec. 

Time  between  exposures  =  2  ysec. 

Shock  Mach  No.  =  12 

p 

Magnetic  Field  Strength  =  1.6  w/m' 


Hall  Parameter 


4.8 
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Figure  3.11(b) 

Image  Converter  Photographs 
of  the  Discharge  Structure 
in  the  Center  of  the  Magnetic  Field  Region 

Exposure  duration  =  0.1  ysec. 

Time  between  exposures  =  1  ysec. 

Shock  Mach  Mo.  =  12 

p 

Magnetic  Field  Strength  =  1.6  w/m' 


Hall  Parameter 


4.8 
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Figure  3.12 


Drawing  of  Frame  #1  of  Figure  3.11(b), 


taken  from  the  original  Polaroid  photograph 
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Figure  3.13 

Angle  between  Perturbation  Current 
and  Original  Current 
vs.  Magnetic  Field 
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3.2.6  Development  of  Turbulence  with  Increasing  Values 
of  Hall  Parameter 

It  has  been  observed7'101  by  workers  using  a  pulsed 
discharge  in  low  temperature  argon  seeded  with  an  alkali 
metal  that,  with  increasing  Hall  parameter,  the  discharge 
becomes  highly  structured  and  very  turbulent.  Workers 
observing  a  discharge  in  a  supersonic  flow  of  low  temp- 
erature  argon  seeded  with  traces  of  sodium  have  reported-3  '  D 
a  similar  increase  in  structure  and  turbulence  with  increas- 
ing  Hall  parameter.  Specifically,  Malikov  has  observed 
that  for  small  oot  two  or  three  broad  pinches  develop.  As 
wt  is  increased,  the  characteristics  become  finer  and  for 
oo t  from  15  to  20  each  pinch  decomposes  into  a  system  of 
finer  filaments. 

A  very  similar  behavior  is  observed  under  the  experi¬ 
mental  conditions  reported  here.  Figure  3.15  shows  photo¬ 
graphs  taken  for  three  values  of  Hall  parameter.  Figure 
3.15(a),  (oot  ^  3),  shows  a  single  somewhat  diffuse  pinch. 
Figure  3.15(b),  (oot  *  4.8),  shows  the  main  pinch  or  streamer 
splitting  up  into  a  number  of  smaller  streamers.  In  add¬ 
ition,  two  secondary  pinches  extend  from  the  upper  electrode 
and  split  into  a  number  of  smaller  streamers.  Figure 
3.15(c),  (oot  =  6.8),  shows  a  great  deal  of  internal 
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structure  to  the  main  streamers.  The  turbulence  inside 
the  main  streamer  causes  the  streamer  to  broaden.  The 
last  frame  of  Figure  3.15(c)  shows  a  general  turbulence 
and  brightening  extending  into  the  background  gas.  Thus, 
it  can  be  seen  that  an  increase  in  Hall  parameter  is 
correlated  with  an  increase  in  structure  and  turbulence 
of  the  discharge. 
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Figure  3.15  (a) 
Variation  of  the  Discharge 
with  Hall  Parameter 

Hall  parameter  =  3 
Exposure  time  =  0.2  y  sec . 
Inter frame  delay  =  2  ysec. 


Structure 


Figure  3.15  (b) 

Variation  of  the  Discharge  Structure 
with  Hall  Parameter 


Hall  parameter  =  4 
Exposure  time  =  0. 
Interframe  delay  = 


.8 

1  ysec. 

1  ysec. 


Figure  3.15  (c) 

Variation  of  the  Discharge  Structure 
with  Hall  Parameter 


Hall  parameter  =  6 
Exposure  time  =  0. 
Interframe  delay  = 


.  3 

2  psec. 

2  ysec. 


CHAPTER  4  BEHAVIOR  OF  THE  CURRENT  STREAMERS 


AFTER  FORMATION 


4.1  Sequence  of  Events  after  Formation 

As  shown  in  Figure  3.11(a),  the  first  noticeable 
feature  is  the  appearance  of  a  number  of  luminous  stria- 
tions  (frames  2  and  3  of  Figure  3.11(a)  )  at  an  angle 
close  to  that  predicted  by  equation  A. 3. 36.  As  the 
striations  develop,  one  or  two  of  them  tend  to  carry  most 
of  the  current  and  these  striations  expand  to  become 
highly  luminous  streamers  oriented  primarily  along  the 
vxB  or  Faraday  direction.  The  secondary  striations  then 
become  less  luminous  and  appear  to  disintegrate  (frames 
4  and  5  of  Figure  3.11(a)  ).  The  primary  streamers  then 
continue  to  increase  in  luminosity  and  break  up  into 
secondary  streamers  at  a  characteristic  angle  to  the 
primary  streamer  as  described  in  Section  3.2.5.  As  the 
current  through  the  primary  streamers  increases,  the 

force  on  the  streamers  also  increases,  decelerating 
them.  This  deceleration  is  dealt  with  in  the  next  section. 
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Further  development  of  the  current  structure  is  limit 
ed  by  the  flow  duration  time  of  the  shock  heated  gas. 

When  the  highly  conductive,  luminous  driver  gas  reaches 
the  electrode  region,  the  vx$  or  Faraday  current  will  be 
conducted  by  the  driver  gas,  thus  removing  the  energy 
source  needed  to  maintain  the  streamers  in  the  shock 
heated  gas.  This  is  clearly  shown  in  Figure  4.1.  The 
driver  gas  (  the  very  luminous  region  at  the  right  hand 
side)  begins  to  intrude  into  the  electrode  region  by  the 
third  frame.  The  process  continues  in  the  fourth  frame 
and  the  streamers  begin  to  decrease  in  luminosity.  By 
the  fifth  frame,  the  streamers  have  almost  completely 
disappeared. 

There  is  some  evidence  to  suggest  that  the  streamers 
are  initiated  from  a  region  of  high  luminosity.  Figure 
4.2  shows  photos  taken  at  right  angles  and  slightly  closer 
to  the  driver  compared  to  Figure  4.1.  The  shock  heated 
gas  flows  from  right  to  left.  The  highly  luminous  region 
on  the  .right  is  the  driver  gas.  The  two  shadows  in  the 
middle  are  the  supports  for  the  magnetic  field  coils. 

An  expanding  luminous  region  is  seen  just  behind  the 
shock  front.  This  luminous  region  is  produced  by  the 
focussing  effect  of  the  applied  magnetic  field.  The 
charged  particles  will  tend  to  follow  the  field  lines 
and  the  gas  will  be  compressed  as  the  field  lines  converge 
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Ref  erring  again  to  Figure  4.1,  it  can  be  seen  that  most 
of  the  streamers  appear  to  emanate  from  this  luminous 
region.  It  can  also  be  seen  that  streamers  emanate  from 
electrode  spots  on  the  upper  electrode  (cathode) .  It 
may  be  that  these  highly  luminous  regions  provide  the 
perturbation  to  trigger  the  instability.  More  investi¬ 
gation  is  required  to  fully  understand  the  role,  of  these 
luminous  regions  in  the  formation  of  the  streamers. 

4.2  Force  Balance  on  the  Current  Streamer  and  the 

Interaction  of  the  Streamer  with  the  Gas  Flow 

In  this  section,  we  will  consider  the  balance  of 
forces  acting  on  the  current  streamer.  The  first  major 
force  is  the  magnetic  interaction  or  jxB  force  which 
acts  in  an  upstream  direction  and  tends  to  decelerate  the 
streamer.  The  second  force  which  must  be  considered  is 
the  hydrodynamic  or  convective  drag  force  which  acts  down¬ 
stream  and  tends  to  increase  the  velocity  of  the  streamer 
towards  the  flow  velocity  of  the  shock  heated  gas.  The 
third  force  is  the  inertia  force  due  to  the  deceleration 
of  the  streamer. 

There  exists  a  substantial  literature  on  electric  arcs 
in  a  transverse  gas  flow  in  the  presence  of  a  magnetic 
field.  References  10,  12,  13,  14,  15,  18,  24,  25,  30,  61, 
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Figure  4 . 1 

Image  Converter  Photographs 
Taken  Through  Center  of 
Magnetic  Field  Coils 


Exposure  duration:  first  frame,  0.5  ysec. , 

second  to  fifth  frames, 0.2 
Inter frame  delay:  frames  1-2,  8  ysec. 

frames  2-3, 3-4, 4-5,  2  ysec. 
B  ~  1.63  v,T/m2;  Mach  No.  =  12 


ysec. 
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Figure  4 . 2 

Photographs  of  Shock  Heated  Gas  as  Shock 
Front  Enters  Magnetic  Field  Region. 

Circular  shadow  at  left  of  last  two  frames  is 
the  shadow  of  the  upper  electrode  (anode) 
Exposure  duration:  0.2  sec. 

Interframe  delay:  2  sec. 

o 

Magnetic  Field  Strength  =  1.63  w/m 


Mach  Number 


12 
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69,  70,  77,  89,  95,  107,  and  118  give  a  representative 

sample  of  some  of  the  more  recent  papers  in  this  area. 

89  lift 

Roman  and  Myers  ,  Winograd  and  Klein  ,  and  Nicolai 
77 

and  Kuethe  ,  among  others,  have  shown  that  the  arc 
behaves  like  a  solid  cylinder  and  is  essentially  impervious 
to  the  gas  flow,  under  their  experimental  conditions. 

Their  experimental  conditions  consisted  essentially  of  an 
arc  subjected  to  a  cross  flow  of  air  or  nitrogen  at  room 
temperature.  The  drag  force  caused  by  the  cross  flow  was 
balanced  by  the  force  produced  by  an  externally 

77  89  118 

applied  transverse  magnetic  field.  The  above  authors  '  ' 

found  that  the  magnetic  field,  B,  required  to  balance  the 
arc  is  approximately  proportional  to  the  square  of  the 
transverse  gas  blowing  velocity,  .  This  dependence  is 
the  same  as  that  obtained  by  equating  the  aerodynamic  drag 
force  and  the  magnetic  retarding  force  on  a  conducting 
cylinder  of  unit  length: 

B  I  =  h  C~  D  p  V2  ...(4.1) 

a  xa  2  D  a  00  00 

where  I  is  the  current  flowing  along  the  axis  of  the 
a 

cylinder,  D  is  the  significant  dimension  (diameter) , 

Cp  is  the  drag  coefficient^8,  and  pm  is  the  free  stream 
gas  density.  The  above  authors77'89'118  found  that  the 
drag  coefficient,  CD,  defined  by  the  above  equation,  is 
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very  close  to  unity  for  a  subsonic  or  near  sonic  cross 
flow.  This  is  approximately  the  value32'35'64  for  a 
solid  cylinder  in  cross  flow  in  the  same  range  of  Rey¬ 
nolds  number.  Figure  4.3  gives  a  graph  of  versus 
Reynolds  number  for  a  solid  cylinder. 

The  value  of  the  drag  coefficient,  CD,  is  a  measure 
of  the  extent  of  interaction  between  the  arc  and  the  flow. 

Other  workers  have  found  higher  values  of  CD  for  a  balanc- 

70 

ed  cross  flow  arc.  For  instance,  Myers  et  al  ,  working 

with  an  argon  arc  in  a  flow  velocity  of  up  to  100  m/sec 

2 

balanced  by  a  magnetic  field  of  up  to  3  w/m  ,  found 
values  of  CD  approximately  six  times  greater  than  that  for 
a  solid  cylinder.  Myers  et  al  correctly  point  out  that 
this  indicates  a  much  more  pronounced  arc-flow  inter¬ 
action  than  in  the  case  of  a  solid  cylinder  of  similar 
size.  Bond  and  Wickersheim  ,  working  with  an  arc  in 
sulphur  hexaf louride,  also  observed  a  more  pronounced 

interaction.  Their  observed  drag  coefficient  again  ranges 

15 

up  to  six  times  that  for  a  solid  cylinder.  The  authors 
observed  a  fine  sawtooth  structure  to  the  arc  column  and 
suggest  that  the  value  of  Cq  and  hence  the  degree  of  inter¬ 
action  are  considerably  affected  by  this  arc  structure. 

The  above  discussion  suggests  that  some  idea  of  the 
degree  of  interaction  between  the  current  streamer  and 
the  shock  heated  gas  flow  may  be  obtained  by  investigating 
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the  magnitude  of  the  drag  coefficient  of  the  streamer. 

Since  the  current  streamer  is  being  decelerated,  we  must 
add  an  inertia  force  term  to  equation  4.1.  Thus, 

Ia  B  ~  55  CD  p»  Da  v~  =  ma  "a  ...(4.2) 

where  I  is  the  streamer  current,  B  is  the  applied  trans¬ 
verse  magnetic  field,  CD  is  the  drag  coefficient,  is  the 
shock  heated  gas  mass  density,  Da  is  the  streamer  diameter, 
V  is  the  relative  flow  velocity  between  the  streamer  and 

oo  u 

the  shock  heated  gas,  ma  is  the  mass  per  unit  length  of  the 
streamer,  and  ua  is  the  measured  deceleration  of  the 
streamer. 

Table  4.1  gives  values  of  the  effective  CD  calculated 
from  equation  4.2  using  the  measured  values  of  the  para¬ 
meters  I  ,  B,  pw,  Da,  V^,  ma,  and  ua  shown  in  the  table. 
Values  of  Reynolds  number  for  flow  past  a  solid  cylinder 
of  diameter  D  ,  and  the  value  of  CD  expected  at  this  Nre 
based  on  solid  cylinder  theory  are  also  given.  pOT  was 
derived  from  the  measured  shock  Mach  number  and  undisturbed 
gas  pressure  and  temperature,  using  the  graphs  given  in 
Appendix  1.  D  ,  V  ,  and  uo  were  measured  from  image 
converter  photographs  taken  of  the  current  streamer. 
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rna  was  assumed  given  by  the  relation 

ma  =  <1'/4)D |  ...  (4.3) 

That  is,  we  assume  the  mass  density  in  the  streamer  to 
be  the  same  as  that  of  the  shock  heated  gas.  The  effect 
of  the  temperature  increase  and  resultant  expansion  will 
be  to  reduce  ma  below  the  value  given  by  equation  4.3. 

Thus,  equation  4.3  specifies  an  upper  bound  for  m^  and 
therefore,  from  equation  4.2,  a  lower  bound  for  C  .  Also, 
the  values  of  maua  given  in  Table  4.1  are  much  less  than 
the  term  IB.  Errors  in  rn  u_  then  have  only  a  small  effect 

d  d  d 

on  the  values  of  Cp. 

The  values  of  effective  CQ  given  in  Table  4.1  show 
that  CD  is  from  two  to  three  orders  of  magnitude  greater 
than  the  expected  value  for  a  solid  cylinder.  This  shows 
that  there  is  a  marked  degree  of  interaction  between  the 
shock  heated  gas  flow  and  the  current  streamer .  The 
possible  reasons  for  this  strong  interaction  will  now  be 
investigated . 
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One  possibility  is  the  existence  of  arc  root  or  elec¬ 
trode  phenomena,  that  is,  the  root  of  the  arc  may  be  tied 
to  a  specific  point  on  the  electrode,  thus  inhibiting  the 
motion  of  the  arc.  Of  course,  the  effect  of  this  electrode 
phenomena  would  be  to  decrease  the  arc  velocity  and  thus 
decrease  the  effective  drag  coefficient.  Thus,  arc  root 
phenomena  will  not  explain  the  observed  large  interaction. 

A  second  possibility  is  that  the  measured  luminosity 
width,  D  ,  of  the  streamer  does  not  correspond  to  the  act- 

Cl 

ual  diameter  of  the  interaction  region.  The  highly  lum¬ 
inous  region  could  correspond  only  to  a  small  core  and  the 
region  of  high  pressure  could  extend  over  a  much  larger 
diameter.  This  would  have  the  effect  of  increasing  the 
hydrodynamic  pressure  proportionate  to  the  increase 

in  D-, .  However,  to  increase  the  hydrodynamic  pressure  by 
the  two  to  three  orders  of  magnitude  required  to  balance 
IB  would  require  that  the  effective  Da  be  larger  than 
the  dimensions  of  the  plasma.  This  possibility  is  then 
ruled  out. 

Another  possibility  is  the  existence  of  a  fine  struc¬ 
ture  to  the  streamers,  on  a  scale  much  smaller  than  that 
disclosed  by  the  image  converter  photographs.  One  can 
conceive  of  the  streamer  being  composed  of  a  large  number 
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of  small  streamers,  each  of  which  is  impervious  to  the  gas 
flow.  The  significance  of  this  is  the  introduction  of 
viscous  effects  to  the  flow  around  these  much  smaller 
impervious  cylinders.  The  Reynolds  number  of  the  flow, 
VJDa/v,  will  be  reduced  proportionate  to  the  reduction  in 
Da,  where  Da  is  the  diameter  of  the  small  impervious 
cylinders,  and  v  is  the  kinematic  viscosity  of  the  shock 
heated  gas.  As  shown  in  Figure  4.3,  the  reduction  in 
Reynolds  number  leads  to  a  large  increase  in  the  drag 
coefficient,  CD,  for  flow  around  a  solid  cylinder,  as  the 
flow  enters  the  viscosity  dominated  regime.  Thus,  if 
the  streamer  fine  structure  was  of  a  small  enough  scale, 
the  drag  coefficient  may  be  large  enough  to  explain  the 
results  of  Table  4.1.  Now,  from  Figure  4.3,  a  drag 
coefficient  of  the  order  of  100  to  1000  implies  a  Reynolds 
number  of  about  0.01  to  0.1.  This  implies  that  the  fine 
structure  of  the  streamers  would  be  of  a  scale  three  to  four 
orders  of  magnitude  smaller  than  the  measured  diameter  of 
the  streamers.  From  Table  4.1,  the  characteristic  length 
of  the  fine  structure  would  be  of  the  order  of  10  ^  to 
10“7  meters.  Now,  this  is  of  the  order  of  the  electron  mean 
free  path  which  implies  that  the  fluid  model  breaks  down 
in  this  length  regime.  Thus,  the  possible  existence  of  a 
fine  structure  to  the  streamers  and  the  introduction  of 
viscous  effects  to  the  flow  do  not  appear  to  explain  the 
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anomalously  large  drag  coefficient. 

Another  possibility  for  explaining  the  large  arc— flow 
interaction  is  the  presence  of  turbulence  in  both  the  shock 
heated  gas  and  in  the  streamer.  Pert^^  has  discussed 
the  production  of  turbulence  in  shock  heated  gases.  He 
found  that  the  transition  to  turbulence  can  be  expressed 
in  terms  of  the  hydrodynamic  Reynolds  number  for  pipe 
flow  based  on  the  properties  of  the  shock  heated  gas. 

He  found  that  under  certain  conditions,  the  energy  of 
the  turbulent  field  can  be  propagated  to  the  shock  front 
and  can  cause  a  disruption  and  instability  in  the  front, 
thus  preventing  the  formation  of  a  region  of  shock 
heated  gas.  Under  the  conditions  of  this  experiment,  the 
Reynolds  number  for  pipe  flow  based  on  the  properties  of 
the  shock  heated  gas  is  about  1500  to  2000,  well  below 
Pert's  critical  Reynolds  number  for  shock  front  instab¬ 
ility  of  9000.  However,  the  value  of  Reynolds  number  for 
this  experiment  is  close  to  the  value  for  the  sustenance 
of  turbulence.  Kaufmann^  gives  this  value  as  about 
2300.  It  is  thus  reasonable  to  expect  that  there  will  be 
some  turbulence  of  the  shock  heated  gas  under  the  conditions 
of  this  experiment,  based  on  hydrodynamic  considerations 
alone.  In  addition,  the  results  of  Section  3.2.6  show  that 
a  large  value  of  Hall  parameter  such  as  exists  in  this 
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experiment  can  induce  additional  large  scale  turbulence. 

Thus,  we  can  reasonably  assume  that  there  exists  a  fairly 

high  degree  of  turbulence  in  the  shock  heated  gas. 

The  problem  of  the  interaction  between  an  electric 

arc  and  a  turbulent  gas  flow  has  been  investigated  by 

•  2  5 

Garosi  and  Bekefi  .  They  made  experimental  studies  of  the 

case  of  an  arc  discharge  struck  in  a  stream  of  flowing  argon 

gas,  the  flow  direction  being  parallel  to  the  axis  of  the 

discharge.  The  Reynolds  number  was  varied  from  zero  to 

6300;  the  flow  then  ranged  from  the  laminar  regime  to 

the  regime  of  fully  developed  pipe  turbulence.  The 
2  5 

authors  found  that  turbulence  in  the  flowing  gas  had 
a  very  profound  effect  on  the  properties  and  behavior  of 
the  discharge.  Among  the  observed  effects  were  high  freq¬ 
uency,  fine  grained  oscillations  of  the  electron  density, 
an  enhanced  loss  of  plasma  particles,  and  an  increase  in 
electrical  resistivity  of  the  arc  plasma.  These  effects 
begin  to  appear  in  the  Reynolds  number  range  of  1000  to 
2000,  similar  to  the  Reynolds  number  range  existing  in 
this  experiment. 

The  most  important  effect  reported  by  Garosi  and 
Bekefi,  from  the  point  of  view  of  application  to  the 
present  experiment,  is  the  existence  of  a  bul^  motion  of 
the  arc,  in  which  the  entire  column  is  displaced  almost 
like  a  rigid  cylinder.  Under  their  experimental  conditions. 
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wit  h  the  flow  along  the  axis  of  the  arc,  the  authors 
reported  that  the  displacement  was  quasiperiodic  with  a 
period  of  the  order  of  one  millisecond.  The  authors 
attribute  this  bulk  motion  to  large  scale  fluctuations 
or  eddies  which  appear  in  the  flow  above  a  Reynolds  num¬ 
ber  of  about  1000. 

It  thus  appears  reasonable,  on  the  basis  of  Garosi 
and  Bekefi's  results,  to  expect  an  enhanced  arc-flow 
interaction  whenever  turbulence  is  present  in  the  flow. 

As  discussed  previously,  there  are  two  mechanisms  acting  to 
produce  turbulence  in  the  present  experiment.  The  first 
is  the  production  of  turbulence  by  hydrodynamic  means 
alone,  whenever  the  Reynolds  number  for  pipe  flow  is  of 
the  order  of  2000  or  greater.  The  second  is  the  mechanism 
discussed  in  Section  3.2.6,  where  a  large  Hall  parameter 
can  produce  large  scale  turbulence.  We  would  then 
expect  a  correlation  between  an  increase  in  arc-flow  inter¬ 
action,  manifested  by  an  increase  in  C^,  and  an  increase 
in  flow  velocity,  U2 /  or  Reynolds  number,  N-^g.  We  would 
also  expect  a  correlation  between  an  increase  in  and 
an  increase  in  Hall  parameter  or  magnetic  field  strength. 

Table  4.2  shows  values  of  drag  coefficient,  CD, 
from  Table  4.1  along  with  corresponding  values  of  Hall 
parameter  and  Reynolds  number,  derived  from  Table  4.1. 
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Table  4 . 2 

Correlation  between 
Effective  Streamer  Drag  Coefficient 
and  Hall  Parameter  and  Reynolds  Number  (pipe  flow) 

of  Shock  Heated'  Gas 


-96- 


It  can  be  seen  that  there  is  a  definite  correlation 
between  increasing  Hall  parameter  and  increasing  drag 
coefficient.  There  is  also  a  slight  correlation  between 
increasing  drag  coefficient  and  increasing  Reynolds 
number,  although  the  narrow  range  of  Reynolds  number 
available  here  makes  it  difficult  to  draw  any  definite 
conclusions . 

In  summary,  it  is  proposed  that  the  large  arc-flow 
interaction  is  due  to  turbulence  induced  in  the  shock 
heated  gas  primarily  by  the  mechanism  discussed  in 
Chapter  3.  The  experimental  data  presented  in  Tables 
4.1  and  4.2  appears  to  support  this  proposition. 

4.3  Helical  Instability  of  the  Current  Streamers 

In  an  attempt  to  maximize  the  resolution  of  the 
electronic  camera  and  obtain  maximum  possible  resolution 
of  the  discharge  structure,  single  photographs  were  taken 
with  the  discharge  imaged  onto  as  large  an  area  of  the 
image  converter  tube  photocathode  as  possible.  The 
resulting  photographs  showed  a  definite  helical  structure 
to  the  discharge.  Figure  4.4  shows  two  typical  photo¬ 
graphs  taken  as  described.  A  "double  helix"  structure 
can  be  seen,  with  both  secondary  filaments  forming 
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Figure  4 . 4 

Photographs  Showing  Helical  Structure 
of  the  Current  Streamers 

Linear  dimensions  are  approximately  0.7 
of  actual  size. 

2 

Mach  Number  =12;  B  =  1.63  w/m  ;  cot  -  5 
The  shock  front  and  background  gas  are  not 
visible  in  the  photo  on  the  right;  the  exposure 
and  development  time  were  varied  to  bring  out 
the  more  intense  portions  of  the  image. 
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helices  about  the  axis  of  the  original  streamer.  In  this 
Section,  the  behavior  and  the  reasons  for  the  formation  of 
this  helical  structure  will  be  investigated. 

One  very  likely  mechanism  for  a  structure  of  this 
kind  arises  from  the  existence  of  a  component  of  magnetic 
field  in  the  direction  of  the  streamer,  parallel  to  the 
axis  of  the  helix.  It  is  easily  seen  that  there  exists 
a  radial  component  of  field,  relative  to  the  axis  of  the 
Helmholtz-pair  field  coils,  due  to  incomplete  flux  link¬ 
age  between  the  two  coils.  This  radial  field  will  then 
provide  a  component  parallel  to  the  axis  of  the  streamer 
discharge,  although  the  magnitude  of  this  field  will  vary 
considerably  over  the  interaction  region. 

An  estimate  of  the  magnitude  of  this  radial  field  can 
be  obtained  in  the  following  way.  We  assume  a  knowledge 
of  the  inductance  of  each  coil  separately  and  of  the 
inductance  of  the  coupled  pair  in  place  beside  the  shock 
tube,  as  shown  in  Chapter  2.  The  inductance  of  a  single 
coil,  L-j_,  was  measured  to  be  5.36x10  ^  Henrys.  The 
inductance  of  the  coupled  pair,  was  measured  to  be 

12 . 2  x10-6  Henrys. 
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Now,  the  total  flux  which  links  both  coils  is84 

$21  +  ^12  =  21  M  ...(4.4) 

where  cf>2]_  is  the  flux  linking  the  2nd  coil  due  to  a  current 
I  in  the  1st  coil,  and  cjjgo  is  the  flux  linking  the  1st 
coil  due  to  a  current  in  the  2nd  coil.  M  is  the  mutual 
inductance  of  the  two  coils  defined  as84 

Ln  =  Lg  +  2M  ...  (4.5) 

Now,  the  total  flux  which  is  produced  by  coil  1  but 
which  diverges  out  of  the  interaction  region  without  link¬ 
ing  coil  2  is  given  by 

$1  “  4*21  =  LXI  -  MI  ...(4.6) 

Similarly,  the  total  flux  produced  by  coil  2  which  does 
not  link  coil  1  is 

<J>2  -  $22  =  L-j-I  —  MI  ...(4.7) 

We  can  consider  the  interaction  region  as  being  bound¬ 
ed  by  a  cd  rcular  cylinder  with  its  axis  coinciding  with  uiie 
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axis  of  the  coils.  The  diameter  of  the  cylinder  is  5  cm., 
(the  inner  diameter  of  the  coils) ,  and  the  length  is 
6  cm.  (the  distance  separating  the  two  coils).  Thus,  the 
field  in  the  direction  parallel  to  the  axis  of  the  coils, 
averaged  over  the  interaction  region,  is 


*21  +  4*12 

ttx  (2.5xl0  ‘2)  2 


. . .  (4.8) 


Similarly,  the  average  field  in  the  radial  direction, 
diverging  from  the  interaction  region,  is 


B 


R 


*2  "  *12  +  *1  ”  *21 

wr— —nr  wy  ■in'  -u 

7rx5xlO-2x6xlC-2 


. . .  (4.9) 


Expressing  the  radial  field  as  a  fraction  of  the 
(measured)  transverse  field,  and  using  equations  (4.4)  to 
(4.9),  we  have 


LM  L1 
3L1  "  % 


0.208  x 


.  .  .  (4.10) 
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Using  the  measured  values  of  and  ,  we  have 


B 


R 


0.366 


.  .  .  (4.11) 


This  value  is  only  an  approximate  average  value  of 
Br.  Br  exhibits  a  considerable  spatial  variation,  in 
particular,  Br  vanishes  on  a  plane  parallel  to  and  equi¬ 
distant  from  the  planes  of  the  two  coils. 

The  mechanism  by  which  an  axial  magnetic  field  can 
lead  to  a  helix  structure  can  be  explained  qualitatively 
as  follows:  as  discussed  in  Chapter  3  and  Appendix  3, 
the  current  filament  breaks  up  into  two  secondary  filaments 
each  at  an  angle  a,  given  by  equation  (A. 3. 36),  to  the 
original.  The  secondary  filaments  then  introduce  a  radial 
component  of  current  with  respect  to  the  axis  of  the 
original  filament.  This  radial  current  interacts  with  the 
axial  or  longitudinal  magnetic  field  to  produce  a 
azimuthal  force  which  acts  until  the  radial  current  vanish¬ 
es,  that  is,  until  the  filament  is  twisted  into  a  helical 
shape  such  that  the  former  radial  component  now  flows 
completely  in  the  azimuthal  direction. 

We  will  now  consider  the  question  of  the  determination 
of  the  radius  and  pitch  angle  of  the  helix  and  the  quest¬ 
ion  of  its  stability.  There  is  a  considerable  literature 
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(references  1,  3,  8,  26,  33,  36,  37,  43,  44,  65,  and  108, 
for  example)  on  the  existence  of  a  helical  instability 
in  the  positive  column  and  the  resulting  enhanced  diff¬ 
usion  due  to  this  instability.  The  theoretical  papers 
(ref.  44,  for  example)  start  from  the  equations  of  cont¬ 
inuity  for  ions  and  electrons  and  the  electron  equation 
of  motion.  A  steady  state  solution  is  obtained  and  this 
solution  is  subjected  to  a  perturbation  of  the  form 
exp  [i  (a)t  +  m0  +  kz)]  where  0  is  the  azimuthal  angle. 

A  critical  magnetic  field  of  the  order  of  a  few  kilogauss 
is  found  for  a  given  mode  (m  =  1) .  In  addition,  the 
frequency  and  wavelength  of  the  helix  at  the  critical 
magnetic  field  are  found^.  However,  this  theory  does 
not  fit  the  experimental  conditions  encountered  here,  for 

the  following  reasons:  first,  the  positive  column  was 

2  3  9  8 

considered  to  be  in  the  Schottky  diffusion  state  '  ; 

that  is,  the  particle  loss  is  due  to  ambipolar  diffusion 
and  the  electron  and  ion  concentrations  are  taken  to  be 
zero  at  the  wail  of  the  discharge  tube.  It  is  also  assum¬ 
ed  that  the  temperature  is  constant  throughout  the  posi¬ 
tive  column.  This  theory  is  only  valid  for  low  pressures 
below  about  5  or  10  torr1.  Above  this  pressure  the  posi¬ 
tive  column  gradually  contracts;  the  thermal  transport 
properties  of  the  gas  and  radiation  loss  determine  the 
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size  of  the  column;  a  radial  variation  of  temperature 

appears  and  increases  the  role  of  recombination  processes 

in  the  loss  of  charged  particles.  In  this  experiment, 

where  the  ionization  process  is  unstable,  the  temperature 

inhomogeneity  is  even  more  pronounced  and,  as  we  have  seen, 

the  current  distribution  is  very  non-uniform.  For  these 

4  4 

reasons,  the  Kadomtsev  instability  theory  is  not  appli¬ 
cable  to  this  experiment. 

To  deduce  the  radius  of  the  helix,  we  make  use  of  the 

idea  that  any  conservative  system  will  tend  to  assume  a 

configuration  of  minimum  energy.  If  we  consider  only  the 

potential  energy  of  the  magnetic  field,  then  this  condition 

97 

is  that  the  integral 


be  a  minimum,  where  B  is  the  sum  of  the  field  and  the 
field  due  to  the  helical  current  filament. 

We  can  approximate  this  situation  by  assuming  B^  to 
be  everywhere  uniform  in  space  and  the  "helix"  to  be 
made  up  of  a  number  of  circular  current  loops.  We 
further  assume  that  the  field  of  each  loop  may  be  considered 
independently  of  the  fields  of  the  other  loops.  Vie  tnus 
consider  a  loop  carrying  a  known  current  immersed  in  a 
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uniform  field  Br/  parallel  to  the  axis  of  the  loop.  The 
problem  then  reduces  to  finding  the  radius  of  the  loop 
which  will  minimize  the  quantity 


W 


Blp> 


2 


dv 


.  .  .  (4.13) 


where  B^p  is  the  field  due  to  a  circular  filament  of 
current,  as  given  by  Jackson41.  The  field  of  such  a 
current  distribution  is  a  very  complicated  function  of 
position,  especially  close  to  the  loop;  thus  the  mini¬ 
mization  of  W  is  a  rather  involved  problem  and  probably  can 
only  be  achieved  by  numerical  techniques.  To  simplify  the 
problem  even  further,  we  note  that  the  field  close  to  the 
current  filament  gives  the  biggest  contribution  to  W 
and  hence  we  need  only  consider  B^p  in  this  region.  As 
the  final  approximation,  instead  of  integrating  B  ever 

O 

all  space,  we  will  minimize  B  at  one  representative  point, 
specifically,  the  point  of  intersection  of  the  axis  of 
the  loop  with  the  plane  of  the  loop.  The  value  of  Blp 
at  this  point  is  given  simply  by 


B 


lp 


2a 


. . .  (4.14) 


where  I  is  the  streamer  current,  and  a  is  the  loop  radius, 
s 
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The  minimum  of  the  quantity  B2  =  (BR  +  Blp) 2  is 

obviously  B-^p  =  -BR.  Thus,  the  characteristic  radius  is 
given  by 


Pol 


oJ-s 


=  B 


2a 


R 


.  .  .  (4.15) 


For  the  conditions  of  Figure  4.4  (Is  -  10 3  amps. 


Br  -  0.366x1.63  =  0.6  w/m2) , 


2a 


2  mm. 


From  Figure  4.4,  the  observed  helix  diameters  are  3  to 
5  mm. ,  so  the  above  rather  gross  approximations  appear  to 
predict  the  helix  diameter  fairly  well. 

r — S 

The  pitch  angle  of  the  helix  is  simply  given  by  the 
angle  a  from  equation  (A. 3. 36),  since  this  is  the  angle 
which  the  secondary  filament  makes  with  the  primary  fila¬ 
ment  or  the  axis  of  the  helix.  From  Chapter  3,  this  angle 
should  be  about  40  to  45  degrees.  From  Figure  4.4,  it  is 
seen  that  the  pitch  angle  of  the  helix  is  quite  close  to 
this  value. 

The  question  of  the  stability  of  the  helix  structure 
will  be  considered  next.  The  structure  is  stable  as  long 
as  Bd  and  I  remain  constant.  However,  if  BR  increases  or 
decreases,  then  the  radius  of  the  helix  will  decrease  or 
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increase  respectively.  If  I  increases,  then  the  radius 
will  also  increase.  In  some  of  the  photographic  obser¬ 
vations  made,  the  variation  of  the  radius  was  observed 
as  Is  increased.  A  pronounced  increase  of  the  radius 
was  observed,  up  to  10  mm.  or  more.  At  this  point,  the 
filaments  appeared  no  longer  able  to  sustain  themselves, 
became  very  diffuse,  and  we re  convected  away  by  the  flow 
of  shock  heated  gas. 


CHAPTER  5  APPLICATIONS  AND  SUGGESTIONS  FOR 


FURTHER  RESEARCH 


5.1  Applications 

The  work  reported  here  has  a  number  of  applications 

to  practical  problems.  The  first  is  to  the  design  of  a 

magnetohydrodynamic  (MHD)  generator  of  the  closed  cycle 

type.  The  proposal^ ' ^  is  to  couple  a  nuclear  reactor 

6  7 

of  the  high-temperature,  gas  (helium)  cooled  type 
to  an  MHD  generator.  However,  due  to  constraints  intro¬ 
duced  by  materials,  safety,  and  economics,  the  maximum 
temperature  of  the  gas  at  the  outlet  from  a  reactor  is 
of  the  order  of  1500  degrees  Kelvin.  This  temperature  is 
insufficient  to  obtain  the  thermal  ionization  ensuring  the 
required  electrical  conductivity  of  the  working  fluid,  even 
with  the  use  of  additives  with  a  small  ionization  potential 
(alkali  metals) . 

To  overcome  this  difficulty,  it  is  necessary  to  create 
a  plasma  in  which  preferential  heating  of  the  electrons 
yj_QlUs  the  required  conductivity.  This  separation  between 
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the  electron  temperature  and  heavy  particle  temperature 
can  be  achieved  by  Joule  heating  of  the  electrons49'51,91. 
The  resulting  nonequilibrium  plasma  is  then  subject  to  the 
ionization  instability  discussed  in  the  preceding  chapters. 
This  instability  has  been  found  to  lead  to  a  turbulent 
plasma  state  with  properties  differing  substantially  from 
those  predicted  by  simple  theory.  Specifically,  the  time- 
averaged  conductivity  has  been  found  to  decrease  by  as  much 
as  an  order  of  magnitude  below  that  predicted  by  uniform 
plasma  theory7'52'101'103.  In  addition,  a  non-uniform 
current  distribution  or  the  restriction  of  current  flow 
to  narrow  filaments  leads  to  the  gas  flow-magnetic  field 
interaction  being  confined  to  only  a  small  portion  of  the 
volume  of  the  MHD  channel,  leading  to  highly  inefficient 
operation . 

It  can  be  seen  that  an  understanding  of  the  onset 
conditions  and  the  behavior  of  this  instability  is  a 
most  important  prerequisite  for  the  successful  operation 
of  a  nonequilibrium  MHD  generator. 

Proposals  have  been  made93  that  the  ionization  instab¬ 
ility  may  be  suppressed  by  operating  the  MHD  plasma  at 
higher  temperatures,  where  Coulomb  collisions  dominate  or 
where  the  additive  is  almost  fully  ionized.  The  signifi¬ 
cance  of  the  work  reported  in  the  preceding  chapters  is 
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that  the  instability  will  not  be  suppressed;  it  will  occur 
in  the  Coulomb  collision  dominated  regime  and  in  the 
absence  of  an  additive. 

These  plasma  stability  problems  do  not  occur  in  com¬ 
bustion  MHD  generators  since  the  required  conductivity  can 
be  attained  by  using  oxygen  enrichment,  air  preheating, 
and  seeding  of  the  working  fluid  (natural  gas,  fuel  oil) 
in  the  combustion  chamber.  Temperature  nonequilibrium 
does  not  occur  to  any  appreciable  extent  due  to  the  much 
more  efficient  energy  transfer  to  the  organic  molecules 
during  electron-heavy  particle  collisions.  Accordingly, 
considerable  success  has  been  achieved  in  the  construct¬ 
ion  and  operation  of  combustion  MHD  generators ,  for 

42  . 

example  ,  the  LORHO  and  AVCO  Mark  V  generators  in  the 

USA  and  the  U-02  installation  in  the  USSR.  Continuous 

operating  times  in  excess  of  50  hours  and  power  outputs  up 

to  43  KW  were  reported  for  the  U-02  installation100. 

The  power  was  delivered  through  a  solid  state  inverter 

system  to  the  Moscow  grid.  The  USSR  is  constructing  a 

pilot  plant  of  about  75  MW  installed  capacity;  25  MW 

is  provided  by  the  MHD  generator.  This  installation, 

known  as  the  U— 25  located  in  Moscow,  was  scheduled  to 

commence  operation  during  1970. 
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Success  has  been  not  nearly  as  impressive  in  the  case 

of  nonequilibrium  closed  cycle  generators,  due  to  the  less 

than  complete  understanding  of  the  processes  which  occur 

in  the  MHD  generator  duct.  Recently,  Decker,  Hoffman,  and 
19 

Kerrebrock  have  described  in  detail  the  operation  of  a 
large,  experimental  nonequilibrium  MHD  generator,  using 
helium  seeded  with  cesium.  Load  currents  at  least  five 
times  those  possible  with  frozen  ionization  were  drawn  but 
these  were  still  far  below  those  predicted  by  uniform  plasma 
theory.  The  authors,  although  assuming  in  their  analysis 
that  a  steady  state  exists  in  the  generator  duct,  state 
this  may  well  not  be  the  case  and  that  the  nonequilibrium 
ionization  may  in  fact  occur  by  the  formation  of  convected 
ionization  instabilities. 

Millionshchikov  (vice-president  of  the  Academy  of 
Sciences,  USSR)  and  his  co-authors0  calculate  the  relative 
costs  of  a  helium-cooled  high-temperature  nuclear  reactor 
coupled  to  an  MHD  generator,  compared  to  the  costs  of  a 
conventional  organic  fueled,  steam  powered  plant,  based 
on  a  power  level  of  1200  MW.  They  find  that  the  specific 
yearly  expenses  are  1.5  to  2.5  times  less  for  the  helium 
cooled  reactor  and  the  cost  of  electric  power,  is  two  to 
three  times  less ,  while  the  capital  investments  are  approx 
imately  equal.  The  authors  suggest  that  the  combination 
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of  a  helium  cooled  reactor  with  both  an  MHD  generator 
and  a  gas  turbine  unit  would  provide  higher  technical  and 
economic  efficiencies  than  any  other  known  power  plant 
system.  Thus,  the  incentives  for  continued  work  on  the 
problems  of  nonequilibrium  MHD  generators  are  great. 

Another  application  of  the  work  of  the  preceding 
chapters  is  to  the  extinguishing  of  a  plasma  arc  such  as 
occurs  in  a  circuit  breaker.  As  shown  in  Chapter  1,  the 
effect  of  the  magnetic  field  is  to  transform  the  current 
distribution  from  a  uniform  current  flow  to  a  current  flow 
constricted  in  narrow  filaments.  An  interesting  question 
is  whether  an  applied  magnetic  field  of  sufficient  magni¬ 
tude  can  have  a  similar  effect  on  a  steady  arc.  Tempera¬ 
tures  in  a  typical  circuit  breaker  arc  are  in  a  range  of 
5  8 

1  to  2  eV  ,  very  close  to  the  temperature  range  in  this 
experiment.  Thus,  the  degree  of  ionization  would  be  approx¬ 
imately  the  same  (5%  to  20%) .  The  non-equipartition  bet¬ 
ween  the  electron  and  ion  energies  in  an  electric  arc  in 
atomic  gases  is  well  established^'^,  and  the  criterion 
Te>Ti  is  satisfied.  The  arc  pressure  would  be  somewhat 
higher  and  the  magnetic  field  strength  would  have  to  be 
increased  accordingly  to  obtain  the  appropriate  Hall 
parameter.  If  similar  phenomena  are  observed  as  were 
observed  here,  then  the  arc  should  become  turbulent  and 
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its  effective  conductivity  should  decrease  considerably 
as  noted  in  references  7,  52,  101,  and  103.  The  simul¬ 
taneous  imposition  of  a  longitudinal  magnetic  field  of 
appropriate  magnitude  may  then  lead  to  a  helical  instab¬ 
ility  of  the  filaments  with  consequent  increasing  of  the 
filament  length  and  disintegration  of  the  filaments  as 
described  in  Chapter  4. 

The  enhanced  gas  flow-arc  interaction  discussed  in 
Chapter  4  also  has  application  to  the  circuit  interruption 
problem.  One  common  method  of  interrupting  an  arc  is  by 

CO  o  o 

the  use  of  a  blast  of  high  pressure  gas  '  .  The  rapidity 

with  which  the  arc  plasma  can  be  removed  from  the  electrode 
region  of  the  breaker  depends  on  the  degree  of  inter¬ 
action  between  the  gas  flow  and  the  plasma,  specifically, 
on  the  magnitude  of  the  effective  drag  coefficient,  CQ. 

The  large  increase  in  the  magnitude  of  CQ  observed  here 
(Section  4.2)  may  well  provide  a  means  whereby  the  ionized 
plasma  may  be  removed  more  quickly  from  the  circuit 
breaker. 

Another  application  is  to  the  development  of  high 

7  8 

power  lasers.  It  has  been  suggested  that  the  plasma 
properties  typical  of  a  nonequilibrium  MHD  generator 
may  be  exploited  to  produce  a  population  inversion  by 
MHD  action  in  a  molecularly  loaded  gas  mixture. 
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7  8 

Nighan  et  al  theoretically  studied  this  possibility  and 
showed  that  conditions  favorable  for  the  production  of  a 
vibrational  population  inversion  in  CC>2  can  be  achieved  in 
a  typical  nonequilibrium  MHD  generator,  for  example, 
using  helium  seeded  with  cesium  and  combined  with  nitrogen 
and  carbon  dioxide.  The  ionization  instability  may  well 
occur  here,  limiting  the  laser  power  obtainable.  A  know¬ 
ledge  of  the  onset  conditions  and  behavior  of  the  instab¬ 
ility  then  becomes  of  considerable  importance. 

5.2  Further  Research 

In  addition  to  the  investigations  into  the  possible 
applications  discussed  in  the  last  section,  there  are 
opportunities  for  further  research  into  the  mechanism  of 
the  ionization  instability,  especially  into  an  improve¬ 
ment  of  the  theoretical  description  of  the  phenomena. 

The  theory  could  be  extended  to  the  non-linear,  large 
amplitude  regime,  perhaps  extending  the  work  of  Gilpin 
and  Zukoski^ / 126 .  The  theory  presented  in  Chapter  3 
assumes  an  infinite  (unbounded)  plasma.  Boundary  effects 
could  be  introduced,  perhaps  along  with  a  consideration 
of  the  large  amplitude  regime,  thus  extending  the  work 
of  Nelson76.  Velikov  et  al115  have  done  a  computer 
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simulation  of  i_he  instability  under  quite  restrictive 
conditions  (Te  constant,  Te>>Tu ,  diffusion  and  transport 
properties  other  than  electric  conductivity  neglected) . 
This  work  could  be  generalized  considerably. 

The  determination  of  the  physical  processes  limiting 
the  size  of  the  current  filaments  is  another  area  for 
investigation.  Rough  calculations  of  convective  heat 
transfer  from  the  filament  to  the  flowing  gas  based  on 
consideration  of  the  filament  as  a  solid  cylinder  provide 
reasonable  agreement  with  measured  Joule  heating  rates. 
However,  assumptions  had  to  be  made  regarding  the  filament 
density  and  temperature.  Further  calculations  could  be 
made  in  this  area,  perhaps  including  the  effects  of  turb¬ 
ulence  on  heat  transfer  rates  and  including  experimental 
data  on  properties  of  the  hot  filament.  Radiation  energy 
losses  may  also  make  an  important  contribution  to  the 

p  Q 

total  loss.  Glushov  et  al  suggest  that  a  solution  to 
the  steady  state  equations  (equations  3.29,  3.31,  and  3.14) 
can  lead  to  a  current  pinch,  and  derive  a  critical  Hall 
parameter  for  the  transition.  However,  the  authors  also 
suggest  that  allowance  for  energy  transport  will  change 
the  critical  condition. 

Various  refinements  can  be  introduced  to  the  theoret¬ 
ical  calculations.  For  example,  the  effects  of  energy 


■ 

V 

«-■ 

,  *  ' 

.  * 

. 


■ 


-115- 


nonequipartition  on  the  electron  distribution  function116 
and  on  the  transport  properties6'17'68  can  be  taken  into 
account.  Also,  the  influence  of  the  Ramsauer  effect  on 
the  temperature  nonequilibrium  can  be  considered20'56 

The  substitution  of  He  for  Ar  as  the  shocked  gas  may 
produce  some  useful  results,  especially  since  helium  is 
proposed  as  the  coolant  for  the  high  temperature,  gas 
cooled  nuclear  reactors  under  consideration  as  a  heat 
source  for  MHD  generators.  The  ionization  potentials  for 
He  are  24  eV  for  Hel  and  54  eV  for  Hell.  Thus,  the  gas 
must  reach  slightly  higher  temperatures  than  in  the  case 
of  argon  for  the  instability  to  appear.  In  addition,  the 
large  energy  difference  between  the  two  ionization  poten¬ 
tials  may  produce  a  temperature  region  in  which  the  instab¬ 
ility  is  damped.  There  may  be  a  temperature  range  center¬ 
ed  about  30,000  degrees  K  where  the  first  level  is  essent¬ 
ially  fully  ionized  but  where  the  second  level  has  not 
started  to  contribute  to  the  ionization  to  any  extent. 

The  large  degree  of  coupling  between  temperature  perturb¬ 
ations  and  electron  density  perturbations  will  then  be 
substantially  removed  and  the  formation  of  the  instability 
may  be  prevented. 

interesting  effect  can  occur  m  an  atomic  gas  with 
only  one  bound  electron,  for  example  hydrogen  or  singly 
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ionized  helium.  The  electron  collisional  energy  loss 
term  will  exhibit  a  maximum  with  increasing  Te.  This 
can  be  seen  by  the  following  argument:  the  electron 
energy  balance  equation  is 

°oEy  =  2(me/ma)  nQ  (3k/2 ) (T0-Tg) v 

For  oQ  =  (nee^ ) / (mev ) ,  and  v  «  ne/Tg/^  (Coulomb  collisions), 
the  energy  input  term,  a0E^ ,  then  increases  as  T0^  v/hile 
the  energy  loss  term  is  proportional  to  n| (Te-Tg) /t|^^ . 

Thus,  when  the  degree  of  ionization  approaches  unity, 
and  ne  saturates,  the  energy  loss  term  reaches  a  maximum 
as  Te  increases  while  the  Joule  heating  term  continues  to 

increase.  Note  that  the  collision  frequency  then  decreases 

_  -j  / p  .  3/2 

as  Te  '  and  thus  the  Hall  parameter  increases  as  Tex  . 

The  effect  is  to  introduce  a  singularity  in  the  steady 
state  solution,  that  is,  above  a  certain  value  of  T  ,  no 
steady  state  solution  exists.  If  the  variation  of  Te  with 
Joule  heating  rate  is  considered,  then  above  a  certain  value 
of  Joule  heating  rate  (or  electric  field)  or,  if  the  Joule 
heating  is  provided  by  relative  motion  of  the  gas  and 
magnetic  field,  above  a  certain  value  of  B,  the  electron 
temperature  will  increase  without  bound  or  until  some  other 
energy  loss  orocess  such  as  radiation  or  thermal  transport 
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begins  to  limit  Te. 

The  exact  physical  meaning  of  this  runaway  condition 
is  not  clear,  that  is,  whether  it  will  tend  to  produce 
more  non-uniformities  in  the  current  distribution  or 
whether  it  will  tend  to  smooth  out  and  stabilize  the 
spatial  distribution.  It  is  interesting  to  consider  the 
behavior  of  the  growth  rate  for  small  perturbations 
(equation  3.39)  in  this  region.  Both  the  terms  wtg 
and  (d  In  Te)/(d  In  ne)  will  increase  although  at  different 
rates,  depending  on  the  state  of  the  gas,  such  as  pressure 
and  magnetic  field.  Calculations  made  in  the  same  way  as 
described  in  Section  3.2.1  except  considering  only  one 
ionization  level  of  the  argon  atom  (equation  3.24  replaced 
by  equation  3.14)  indicate  that  the  growth  rate  increases 
very  sharply  in  the  region  of  the  ne  saturation.  This  may 

4  9 

indicate  that,  contrary  to  proposals  made  m  the  literature  , 

the  instability  may  not  be  suppressed  as  complete  ioniz- 

,  .  4  7,101,102,117 

ation  is  approached.  Experimental  results 

made  in  electron-atom  collision  dominated  plasmas  indicate 
that  stability  may  be  increased  as  the  degree  of  ionization 
of  the  additive  approaches  unity.  However,  this  may  not 
the  case  for  a  Coulomb  collision  dominated  plasma. 
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In  conclusion,  the  ionization  instability  dealt 
with  in  this  work  is  a  basic  instability  of  a  partially 
ionized  plasma  in  crossed  electric  and  magnetic  fields 
and  can  manifest  itself  in  MHD  generators,  electric  arcs, 
lasers  utilizing  a  gas  flow  and  a  magnetic  field  for  arc 
balancing,  MHD  lasers,  and  magnetoplasmadynamic  accelerators, 
for  example.  Thus,  continuing  investigations  into  the 
properties  of  this  instability  are  certainly  justified. 
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APPENDIX  1  EQUILIBRIUM  PROPERTIES  OF  SHOCK 

HEATED  ARGON 


Figures  A. 1.1  to  A. 1.3  show  density  ratio,  flow 
velocity,  temperature,  and  degree  of  ionization  of 
shock  heated  argon  in  thermal  equilibrium.  The  curves 
were  derived  from  the  tables  of  de  Leeuw  .  The  initial 
(unshocked)  gas  was  assumed  at  a  temperature  of  300  degrees 
Kelvin  and  a  pressure  of  0.5  mm.  Hg.  is  the  ratio  of 

the  shock  speed  to  the  speed  of  sound,  aj_,  in  the  un¬ 
shocked  gas  (a1  =  3.3*102  m/sec).  P2/Pp  is  the  ratio  of 
the  shocked  gas  density  to  the  unshocked  gas  density. 
u2/(aiMi)  is  the  ratio  of  the  shocked  gas  flow  velocity 
to  the  shock  wave  velocity,  both  velocities  measured  in 
the  laboratory  frame.  T 2  and  a  2  are  the  temperature  and 
degree -of  ionization  respectively  of  the  shock  heated  gas. 


r- 


‘ 

. 


;.'cow»t»tm  xMAOtteK.  wntsttiitonswa&XM 


12  14  16 


18 


Mac h  N umb e  r 


Fiqure  A. 1.1 


Density  Ratio  and  Flow  Velocity  vs.  Mach  No 


1.3 


‘igure  A.  1.2 


Eaui librium  Temperature  of  Shocked  Gas 


vs . 


Mach  Number 


M1 

Figure  A. 1.3 
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APPENDIX  2  RAIZER'S  METHOD  FOR  CALCULATION  OF  AN 

APPROXIMATE  IONIC  CHARGE  IN  THE  REGION  OF 
MULTIPLE  IONIZATION. 


An  approximate  method  of  calculation  of  the  ionic 
charge  has  been  developed  by  Raizer^^.  The  method  rests 
on  two  assumptions.  The  first  assumption  is  that  the 
ion  number  density  nm  and  the  ionization  potentials  I  +1 
are  considered  to  be  continuous  functions  of  the  ionic  charge 
multiplicity  m,  obtained  by  connecting  the  discrete  values  of 
n  and  Im+^  by  continuous  curves.  We  denote  these  two 
continuous  functions  by  n (m)  and  I (m)  respectively.  The 
system  of  recurrence  relations  defined  by  the  Saha  equat¬ 
ion  can  then  be  transformed  into  a  differential  equation 
for  the  function  n (m)  by  replacing  the  finite  differences 
by  differentials: 

n  (m  +  1 )  =  n (m)  +  dn/ dm  ...  ( A . 2 . 1 ) 


Then,  from  Saha's  equation, 
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. . . (A.2.2) 


-146- 


where 


A  = 


27rm_k  3/2 

V3 


and  using ,  from  equation  (A. 2.1), 


n  ( ml- 1 )  d  In  n(m) 

-  =  1  +  — - ■— 


n  (m) 


dm 


1  + 


d  In  n(m) 


dm 


n. 


A 


,3/2 


exp 


I (m+1) 


kT 


. . . (A. 2 . 3) 


We  assume  the  ratio  of  the  partition  functions  is 
approximately  equal  to  unity.  For  argon,  (um+l)/um 
6  for  m=0  ,  1.5  for  m-1,  0.44  for  m=2,  and  2.25  for  m=3 
(see  Table  A. 2.1).  However,  the  variation  in  this  term 
is  small  compared  to  the  variation  in  the  exponential  term. 

The  conservation  of  the  number  of  ions  £nm  =  ng 
becomes 


. . . (A. 2.4) 


where 


is  the  total  number  density  of  ions  and  electrons. 
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The  conservation  of  charge  Ym  nm  =  n^  becomes 

^  III  6 


J m  n  (m)  dm 


n 


e 


. . . (A.2.5) 


where  ne  is  the  total  number  density  of  electrons. 

So  the  average  number  of  free  electrons  per  original 
atom  or  the  average  heavy  particle  charge  multiplicity  is 


. . . (A.2.6) 


The  second  approximation  is  that  the  average  value  of 
the  ionic  charge  multiplicity,  m,  is  equal  to  that  value  of 
m  for  which  the  distribution  function  n  (m)  has  -*a  maximum. 
Obviously,  this  assumption  is  more  justified  the  sharper 
and  narrower  is  the  peak  of  the  distribution  n (m) . 

(Raizer  shows  that  the  peak  is  in  fact  narrow  and  sharp) . 

Thus,  from  equations  (A. 2, 3) and  (A.2.6) ,  we  have  the 
following  simple  transcendental  equation  in  m,  noting  that 
(dn)/(dm)  =  0  at  the  peak  or  maximum: 


I 


.  .  .  (A.  2 . 7) 


where  I  is  the  ionization  potential  of  ions  with  an  average 
charge  multiplicity  m. 
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A  certain  degree  o  f  arbitrariness  in  defining  I 
occurs  which  is  connected  purely  with  the  formalism 
(in  the  exact  theory)  of  assigning  subscripts  to  ioniz¬ 
ation  potentials.  If  we  denote  the  ionization  potentials 
of  an  m-ion  by  Im+^  (the  ionization  potential  of  a  neutral 
atom  is  1^) ,  then  formally  we  should  have  set  1=1 (m+1). 
Sometimes,  however,  the  ionization  potential  of  an  m-ion 
is  denoted  by  I  (the  ionization  potential  of  a  neutral 
atom  is  then  Iq) •  In  this  case,  in  the  Saha  equation,  Im+p 
should  be  replaced  by  lm  and  I  should  be  formally  set  equal 
to  I  (m)  . 

Of  course,  if  we  consider  heavy  elements  at  high  temp¬ 
eratures,  where  the  degree  of  ionization  is  very  high,  the 
arbitrariness  does  not  result  in  any  substantial  change  in 
the  value  of  in  (since  in  this  case,  I  ^  -  Im  <<Im)  . 

Where  the  average  charge  of  the  ions  is  not  large,  however, 
the  arbitrariness  has  a  marked  effect  on  both  values  of 
m  and  values  of  the  thermodynamic  functions.  This  results 
from  the  approximation  of  replacing  discrete  values  by 
continuous  functions. 

Comparison  of  the  results  of  approximate  and  exact 
calculations  shows  that  best  agreement  is  obtained  when  we 
denote  the  ionization  potential  of  an  m-ion  by  Jm+1 ' 
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setting  Iq  =  1(0)  =  0,  but  referring  the  "average"  value 
of  the  potential,  I,  to  the  point  in+^,  i.e.,  assuming 
I  =  I (m+%) .  This  procedure  appears  to  be  quite  natural, 
when  we  note  that  the  sequence  of  discrete  values  of  m  is 
separated  by  finite  intervals  Am  =  1. 

Now,  for  argon,  the  function  I (m)  can  be  approximated 
very  closely  (see  Table  A. 2.1)  for  the  first  four  or  five 
stages  of  ionization  by  the  expression 

I (m)  =  15  m  electron-volts  ...(A. 2. 8) 


Taking  the  logarithm  of  equation  (A. 2. 7),  we  have 


I  (in+'S) 


k  Te  In 


AT 


3/2 


mn 


g  _ 


. . . (A. 2.9) 


which  can  be  easily  solved  numerically  to  obtain  m  and 
thus  ne  for  a  given  Te  and  n  . 
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Table  A.  2.1 

Table  of  Ionization  Potentials 
and  Partition  Functions  (to  first  order) 
for  Argon  (from  ref.  2) 


Ionization 

potential 

(eV) 

Partition 

Function 

(v 

AI  (atom) 

15.755 

1 

All  (singly- 
ionized) 

27.62 

6 

All  I 

40.90 

9 

AIV 

59.79 

4 

AV 

75.0 

9 

AVI 

91.3 

6 

AVI  I 

124.0 

1 

AVI  1 1 

143.46 

2 

AIX 

422.6 

1 

AX 

479.4 

6 

AX  I 

538.9 

9 

APPENDIX  3  DETAILED  CALCULATION  OF  THE  GROWTH  RATE 


EXPRESSION  FOR  IONIZATION  INSTABILITY 


The  following  derivation  is  an  extension  of  the  work 


Their  work  is  extended  to  the  multiply-ionized.  Coulomb 
collision  dominated  regime. 

The  procedure  is  to  consider  small  perturbation 
solutions  to  the  three  equations : 

i)  Generalized  Ohm's  Law  (eq.  3.21) 


a 


0 


ii)  Electron  energy  balance  equation  (eq.  3.23) 


dne 

MPXS'lMi 


dt 


j • E  -  2(me/ma)nev(3k/2) (Te-Tg)  ...(A. 3. 2) 


iii)  Ionization  equilibrium  expression  (eq.  3.25) 
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. . . (A. 3.3) 
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We  assume  small  initial  perturbations  to  the  equil¬ 
ibrium  values  of  n0,  T0,  t,  and  j  as  in  equation  (3.26), 

We  take  the  mean  free  time  for  electron-ion  collisions  (eq.3.7) 

as  i 

32/77  (kTe)3/2 

T  =  - - - - — — — - —  ...(A.  3. 4) 

e4  Z  ne  InA 

3  /2 

Thus,  using  t  «  T0  ~/ne  (neglecting  the  InA  variation), 
and  inserting  the  perturbed  expressions  for  nQ  and  T0 ,  we 
have 


T  =  T 


3T 

1  +  — 
2T 


el 


eO 
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el 

1*3  sown 

ne0 


. . . (A. 3 . 5) 


Similarly , 


v  =  1/t  = 


1/ t 
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1  + 


nel  3Tel 


ne0  2TeO 


... (A. 3 . 6 ) 


2  3/2 

Using  a  =  (n  e  )/(mev)  or  a  Te  ,  we  have 


a  =  a. 


3T 


1  + 


el 


2T 


eO-* 


. . . (A. 3.7) 


Inserting  the  perturbed  expressions  for  n0,  Tq,  E, 


and  3",  Ohm's  Law  becomes,  to  zero  th  order. 
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o  E 
o  o 


*o  +  “T03ox(B/B) 


.  .  .  (A. 3. 8) 


We  neglect  the  diffusion  terms  in  Ohm's  Law, 
To  first  order.  Ohm's  Law  is: 
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(A. 3. 9) 


The  energy  balance  equation  becomes,  to  zero'th 
order , 


Jo  o 


=  (2me/ma) ( 3k/2 ) (ne0/TQ) (TeQ-T  ) 


(A. 3. 10) 


To  first  order,  the  energy  balance  equation  is 
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Note  that  including  the  diffusion  terms  in  Ohm's 
Law  would  add  imaginary  terms  to  Ohm's  Law  (equation  A. 3. 9) 
and  thus  on  the  right  side  of  the  dispersion  relation 
(A. 3. 11),  through  the  Joule  heating  terms  jg’EQ  and 
These  terms  would  not  affect  the  growth  rate  but  would 
introduce  a  finite  phase  velocity  (see  refs.  34  and  72) . 

The  ionization  equilibrium  expression  (A. 3. 3)  becomes, 
by  inspection. 
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We  note  that,  from  equation  (A. 3. 3), 
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We  need  two  more  conditions  in  order  to  eliminate  the 
first  order  quantities  from  the  dispersion  relation  (A. 3. 11). 
We  shall  exclude  electromagnetic  wave  propagation  from 
consideration  and  assume  that  the  electric  field  has  a 
potential  p  so  that  $  =  -gradcf  and  thus 


curl  Is  =  0 


. . . (A. 3. 13) 


We  also  make 


the  plasma  approximation  that  the  net  charge 
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density  is  zero  and  thus 
div  1  =  0 


. . . (A. 3.14) 


From  div  j  -  0  and  the  form  of  j,  3^’K  =  0.  Thus, 
\j±x^\  =  jpK. 

From  curl  E  =  0  and  the  assumed  form  of  E,  E^xK  =  0, 
and  thus  2-^K  =  E^K. 

Thus,  since  j^-LK  and  E-jMK,  must  be  -L  to  2^  or 

3i,Ei  =  °- 

So,  we  have 


-  0  E^K 

,t  ^  i  .  ->  -> 

blxK|  -  DXK  Ex-K 


0 

EiK 


■f 
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•  ...  (A. 3 . 15) 


We  see  that  there  are  two  possible  orientations  of  the 


and  K  vectors : 


Figure  A. 3.1  Figure  A. 3. 2 

( tt/4  )  <(p<  (tt/2) 
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Figure  A. 3.1  corresponds  to  the  perturbed  current 

vector  3^  on  the  electric  field  side  of  the  unperturbed 

current  3^.  Figure  A. 3.2  corresponds  to  3^  on  the  side 

°f  Dq  opposite  to  <J>  is  the  angle  between  3^  and 

6  is  the  angle  between  j.  and  Ert. 

0  0 

We  will  seek  to  express  all  products  of  vectors  in 
terms  of  or  iox^‘  In  the  following  expressions,  when 

a  double  sign  is  written  (±  or  J),  the  upper  sign  refers 
to  Figure  A. 3.1  and  the  lower  to  Figure  A. 3. 2. 

Next,  we  evaluate  the  terms  5q*^i  an<3  5]_*Eq  ec3uatl°n 
(A. 3. 11).  From  equation  (A. 3. 15) and  Figures  A. 3.1  and  A. 3. 2, 

E 

30«21  =  (30«S)  — =  j  qe1  cos  $  ...(A. 3. 16) 

K 

Jl-So  =  -t  I  KxEq  I  —  for  Figure  A.  3.1 

K 

Now,  for  Figure  A. 3. 2,  Kx£Q  changes  sign  as  ux  increases. 
}rfio  goes  from  positive  to  negative  as  wx  increases. 

Now,  from  the  zero-order  Ohm's  Law,  equation  (A. 3. 8), 

a  E  *K  =  t,**  +  - tf0*B)*K 

o  0  u  g 

and  (}Qxg)x£  =  (]q-K)S 

perpendicular  to  S. 


since  all  quantities  are 
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So , 
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From  Figure  A. 3.1, 
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From  Figure  A. 3. 2, 
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Thus,  jg’Eo  is  given  by 
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(sin  $  ±  cotq  cos  )  ...(A. 3. 18) 


We  will  need  Eq*K  in  terms  of  K  and  We  again  use 
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So, 
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IS* Jo  I 


(OX 


jQ  K  cos  <p  “ 


j  qK  sin  <f> 


. . . (A. 3. 19) 


Next,  we  seek  to  express  5^  and  15^  in  the  above 
expressions  in  terms  of  ne^,  Te^,  <j> ,  and  the  zero-order 
quantities . 

To  find  we  take  the  outer  product  of  ^  with  the 

first-order  Ohm's  Law,  equation  (A. 3. 9): 

Term  1  vanishes  from  equation  (A. 3. 15). 

Term  2  is,  from  Figures  A. 3.1,  A. 3. 2,  and  equation  (A. 3. 17), 
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3  Te^  3  Te^ 
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Term  3  is,  from  equation  (A. 3. 15), 


=  +j1K(S/B) 

Term  4  vanishes  since  (^xS)  XK  =  (3]_*K)S  -  (S •  S )  □  j_ ,  and 

=  0  from  equation  (A.  3. 15)  and  B*K  =  0  since  ^  is 

perpendicular  to  all  other  quantities. 
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Term  5  is 


(jQxB)xK  =  (30-K)B  -  (B-K)j0, 


and  S-S  =  0  by  the  s 
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ame  reasoning  as  in  term  4 .  So  term  5  is 
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l2T 


eO  n 


OqKCcos  <j)  )  (§/B)  ,  and  the  first- 


e0> 


order  Ohm's  Law  becomes 


or 


3T 


el 


2T 


tjQK  sin  <J>  +  U3T0joK  cos  ^ 


eO 

=  +  a)T, 


] 


3T^t  n 


2T 


el  _  A1el 
eO  ne0 


jQK  cos  (p 


Ji  = 


3T 


L2Te0 


el. 

— j  sm  (p 


0 


+ 


n 


el 
LneO 


GOT 


o  00  COS  * 


. . (A. 3 .20) 


To  find  ,  we  take  the  inner  product  of  the  first  - 
order  Ohm's  Law,  equation  (A. 3. 9),  with 


3Tel  "To 

°oElK  +  °o“r"  V*  = 


2T 


eO 


B 


(3xxg)  •£ 


0)  T 


B 


3Tel  nel 

~y  ~y  -y 

(J0xB) -K 

-3TeO  ne0- 

- 

Now,  from  Figures  A. 3.1  and  A. 3. 2  and  equation  (A. 3. 15) , 


$•(]]_»§)  =  S-  *  0,  and 


->  ~y . 

K*  ( D  ]_ XB ) 


=  BKjx 
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£•  (D0xS)  =  S-  (S><30)  *  0 

(?oxS>  I  =  lsMM0l 

So,  the  first  order  Ohm's  Law  becomes,  using  equation  (A. 3. 19), 


3T 


°oEl  + 


el 


2T, 


b  (‘ 


e0 

;03Tojl  +  03T 


3q  (cos  <}>  + 

3T 


utq  sin  <J> 


) 


el  nel\  . 


o 


2T 


eO  ne0, 


j0  sin  <j> 


or 


3Tel  .  A  .  nel  . 

aQE-^  +  — —  3  q  cos  <f>  _j.  cot 


2T 


j0  sin  <J> 


eO 


n 


eO 


-“Vi  =  0 


Using  equation  (A. 3. 20)  for  3^, 


o  l 


n. 


3Tel  .  3Tel 

+  — —  j  Q  cos  <j>  _  — ■ —  wtq3  0  sin  <f> 


2T 


eO 


2T 


eO 


el  2  .  . 

+  *" - ~  (  CO  T ©  )  jQ  COS  4> 

ne0 


n 


el 


COT 


jQ  sin  <j)  =  0 


n 


eO 


Simplifying , 


£0  3Tel 
ao  2TeO 


n 


|-cos  <j>  ”  to t 0  sin 


+  —  co t o  —  (-0)To  COS  <j>  t  sin  4) 

ao  °  ne0 


) 


. . . (A. 3.21) 
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i  :  ■ 
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Combining  equations  (A. 3. 18)  and  (A. 3. 20), 
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i  2 
30 


u0 
nel 


3Te1  .2  nel  o  9 

■  '  “  sm  (j>  +  — ■ —  (03T  )  COS  ^  <p 

2TeO  ne0  ° 
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eO  n 


e0> 


. . . (A. 3. 22) 


Combining  equations  (A. 3.16)  and  (A. 3. 21), 


Vfii  = 


-i  2 

3o 
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3T 


2T 


el  o 
cos  <p 

eO 
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el 


2  2 

(ootq)  COSZ4> 


n 


eO 


3T 


el  nel 


+ 


2T 


0)to  cos  tj)  sin  <j> 
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. . . (A. 3 .23) 


(A. 3. 22) 

and 

(A.  3. 

-f  r> 

3 1  *  E0  = 

ii 

3Tel 

ao 

2TeO 

£ 


+  n  nel 

2  — —  cot  sin  (h  cos  <f> 
o 


sin  (j>  - 


COS^I 


n 


eO 


Thus,  the  dispersion  relation,  equation  (A. 3. 11),  becomes. 


lfflne0kTe0 


3n 


el 


3T 


el 


dl  n 


+ 


el 


n 


+  I 


el 


2n 


eO 


2T 


eO 


ngdm  kTe0 


ne0kTe0 


i2 
3  0 


o 


3T 


el /  .  2, 

[  sm  <p 


2TeOV 


2  \  i  nel 

cos  <p)  _  2 - cos  <J>  sm  cp 


n 


mm  nkT  „ 
e  eO  eO 


m„  x 
a  o 


n 


el 


eO 


\  /  T  \  T  ,  nQl  /  T 

e1'  i  -  -M+-£i  +  —  fi  -  -a. 


ne0  2TeO 


Teo/  Te0  ne0 


eO 


-3 


♦ 


*  " 


. 


-162- 


Now,  from  equation  (A. 3. 10), 


me  kTe0 


nu  x 
Cl  o 


J0  ^0 
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1  - 
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eO 
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eO 


From  equation  (A. 3. 8),  jQ  e( 


io/°, 


Thus 


me  kTe0 
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aone0 
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eO 


Also,  from  equation  (A. 3. 12), 


el 


eO 


d  ln  Te\nel 

trmw*  —  —  ■■  I 

‘d  ln  ne/ne 0 


The  dispersion  relation  becomes 
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320 
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aone0I  (2  d  ln  ne 


sin^(j>  -  cos^<}>]  l2ooT0sm 


—  2  + 


3  d  ln  T 


2  d  In  n 


Te0  \d  ln  Te 


,T  n-T  /d  In  n 
eO  g'  e 


3  kT 


eO 


[1  + 


d  ln  T  \  n  ~  1  dl 

e  \  +  eO  _ _ +  ^ 


-1 


2  I 


d  ln  n. 


n_  I  dm 

g 


<{>cos<J> 


(A. 3. 24) 


We  see  that  the  amplitude  of  the  growth  rate  depends 
on  (f> ,  the  angle  between  and  IL  This  means  that  there 
is  a  preferred  directional  dependence  of  the  instability. 


.  •  *  • 
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To  find  the  maxiraum  value  of  the  growth  rate,  iw, 
and  the  angle,  <f> ,  for  which  the  growth  rate  is  a  maximum, 
we  maximize  the  function 

3  d  In  Te  /  \ 

F  =  «_  — —  (sin^cj)  -  cos2<}>]  t2o)i0  sin<|)Cos<}> 

2  d  in  ne  \  / 


with  respect  to  <|> .  That  is,  we  solve  for  <$>max  and  Fmax 
under  the  conditions  ( dF )  /  ( d <J> )  =  0;  (d2F)  /  (dcf>2 )  <  0. 

Now, 

3  d  In  Te 

p  = —  — - — —cos  2  4>  ~  0)Tosin2({)  ...(A.  3. 25) 

2  d  In  nQ 


dF  d  In  T 

— —  =  3  »— — — — ~  sin2t})  2u3t0cos2(|)  .  .  .  (A.  3 . 26) 

d<+>  d  In  ne 

d2F  d  In  T 

_ _  6_ - - —  cos2<j>  7  4o)t  sin2c{>  ...(A. 3. 27) 

O  T  i-  O 

d <J>  d  In  ne 

The  condition  (dF)/  (d<j))  =  0  leads  to,  from  equation  (A.  3. 26)  , 


tan 

2*  =  ; 

WTo 

3 

d  In  Te 

.  2 

d  In  nQ  _ 

The  condition 

(d2F)/(dcj>2) 

<  0  leads 

. . . (A. 3. 28) 


sin  2<f  >  0;  cos  2*  <0.  Thus,  U/2)  <  24,  <  it. 
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Thus , 


cos  2  d> 

max 


and , 


• • • (A. 3 . 29) 


. . . (A. 3. 30) 

and  Fmax  is  given  by,  from  equations  (A. 3. 25),  (A. 3. 29), 
and  (A. 3.30)  , 


9  2 

Similarly,  for  Figure  (A. 3. 2),  the  condition  (d  F)  /  (d<J>  )  <  0 
leads  to,  from  equation  (A. 3. 27),  cos  2<j>  <  0;  sin  2<j>  <  0. 
Thus,  it  <  2$  <  ( 3 tt/2 )  ,  and 


cos 


26 

^max 


3  d  In  T 
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2  d  In  n 

e 


(wtq)2  + 


3  d  In  T 


2  d  In  n 


.  . .  (A. 3. 32) 


-165- 


sin 


26 

max 


. . . (A. 3. 33) 


and  Fmax  is  given  by,  from  equations  (A. 3. 25),  (A. 3. 32), 
and  (A. 3 . 33)  , 


Thus,  we  see  that  F  is  the  same  for  both  orient- 

max 

•  ^  — v 

ations  of  3^  with  respect  to  jQ,  given  by  Figures  A. 3.1 

(equation  A. 3. 31)  and  Figure  A. 3.2  (equation  A. 3. 34). 

This  means  that  the  growth  rate  is  a  maximum  for  both 

->  28 

orientations  of  3^  with  respect  to  3q  and  implies  that 
a  streamer  will  break  up  into  two  parts,  each  at  the  charac¬ 
teristic  angle  to  the  original. 

The  maximum  value  of  the  growth  rate  is : 
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For  both  orientations  of  j ^ , 


2  4>  =  -  tan 


-1 


of  j 

1' 

0)  T 

O 

- 

3 

d  In 

T 

e 

2 

d  In 

n 

where  the  value  of  tan  1  is  taken  between  0  and  u/ 2. 

Thus,  the  angle  between  jQ,  the  average  or  unperturbed 
-> 

current,  and  j-^,  the  perturbation  current,  is  (see  Figures 
A.  3 . 1  and  A. 3 . 2)  : 


IT 
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h  tan 


-1 


WT 


3  d  In  T 


2  d  In  n 


e  -* 


.  .  .  ( A .  3 . 3  6 ) 


where,  again,  tan  ^  is  taken  between  0  and  tt/2 


